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CUBATURE FORMULAS, GEOMETRICAL DESIGNS, 
REPRODUCING KERNELS, AND MARKOV OPERATORS 

O ' 

o 

d '. Pierre de la Harpe and Claude Pache 

Abstract. Cubature formulas and geometrical designs are described in terms of reproducing 
kernels for Hilbert spaces of functions on the one hand, and Markov operators associated to 
orthogonal group representations on the other hand. In this way, several known results 
for spheres in Euclidean spaces, involving cubature formulas for polynomial functions and 
spherical designs, are shown to generalize to large classes of finite measure spaces (fl,cr) and 
appropriate spaces of functions inside L 2 {VL,a). The last section points out how spherical 
designs are related to a class of reflection groups which are (in general dense) subgroups of 
orthogonal groups. 
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Let (fi, a) be a finite measure space and let T be a vector space of integrable real- valued 
functions on O. It is a natural question to ask when and how integrals J n (p(uj)da(uj) can 
be computed, or approximated, by sums J2 xeX W(x)(p(x), where X is a subset of O and 
W : X — > R* + a weight function, for all ip G T. 

When Q is an interval of the real line, this is a basic problem of numerical integra- 
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tion with a glorious list of contributors: Newton (1671), Cotes (1711), Simpson (1743), 
Gauss (1814), Chebyshev (1874), Christoffel (1877), and Bernstein (1937), to quote but a 
few; see [Gauts-97] for some historical notes, [Gauts-76], and [DavRa-84]. The theory is 
inseparable of that of orthogonal polynomials [Szego-39]. 

When Q is a space of larger dimension, the problems involved are often of geometrical 
and combinatorial interest. One important case is that of spheres in Euclidean spaces, 
^ with rotation-invariant probability measure, as in [DeGoS-77] and [GoeSe-79]. Work 

related to integrations domains with dim(fi) > 1 goes back to [Maxwe-77] and [Appel-90]; 
see also the result of Voronoi (1908) recalled in Item 1.15 below. Examples which have 
been considered include various domains in Euclidean spaces (hypercubes, simplices, 
Item 1.20), surfaces (e.g. tori), and Euclidean spaces with Gaussian measures (Item 3.11). 
There are also interesting cases where O itself is a finite set [Delsa-78] . 

Section 1 collects the relevant definitions for the general case (O, a). It reviews several 
known examples on intervals and spheres. 

Our main point is to show that there are two notions which are convenient for the 
study of cubature formulas, even if they are rarely explicitely introduced in papers of this 
subject. 

First, we introduce in Section 2 the formalism of reproducing kernel Hilbert spaces 
([Arons-50], [Krein-63], [BeHaG-03]), which is appropriate for generalizing to other spaces 
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various results which are standard for spheres. See for example the existence theorem for 
cubature formulas (X,W) with bound on \X\ (Items 2.6-2.8), Propositions 2.10 & 3.4 
on tight cubature formulas being geometrical designs, Proposition 2.12 on the set of dis- 
tances Dx = {c £ ]0, oo[ | c = d{x, y) for some x, y G X, x ^ y} if O is a metric space and 
if "Condition (M)" is satisfied, and Items 2.12-2.13 for a single equality which guarantees 
that a pair (X, W) is a cubature formula of strength k. 

Section 3 concentrates on spheres, on consequences of the existence of the central sym- 
metry, and in particular on spherical designs of odd strengths. An example of "lattice 
construction" (Item 3.10) provides a cubature formula in S 7 of strength 11 and size 2400. 
Item 3.11 hints at some very rudimentary constructions of Gaussian designs. 

Secondly, we consider in Section 4 Markov operators (called elsewhere "convolution 
operators", "difference operators", "Hecke operators", or "discrete Laplace operators"). 
On the one hand, they provide an alternative definition of cubature formulas on spheres 
[Pache-04], and at least part of this can be generalized to other Riemannian symme- 
tric spaces of the compact type (work in progress); also, they show a natural connection 
with the work of Lubotzky, Phillips, and Sarnak (see Theorem 4.7). On the other hand, 
they suggest an interesting class of examples of infinite reflection groups, as shown in 
Section 5 where we state and prove an unpublished result of B. Venkov on remarkable sets 
of generators of the orthogonal group 0(8, Z[l/2]). 

The subject of cubature formulas has natural connections with several other subjects. 
We give brief mentions to: 

o representations of finite groups (see [Banna79], [HarPa-04], and Example 1.14), 
o lattices in Euclidean spaces (see [Venko-84], [VenMa-01], and Examples 1.16 & 
3-10), 

o Lehmer's conjecture (Example 1.16) and Waring's problem (Item 3.9) from number 
theory, 

o Dvoretzky's theorem from Banach space theory (see [LyuVa-93] and Proposi- 
tion 3.12). 

It is a pleasure to thank Boris Venkov, as well as Bela Bajnok, Eiichi Bannai, Gael 
Collinet, and Martin Gander for many conversations and indications which have been 
most helpful during this work. 

1. Cubature formulas, designs, and polynomial spaces 

Let (O, a) be a finite measure space. By a space of functions in L 2 (0, cr), we mean a 
linear subspace of genuine 1 real- valued functions on O, given with an embedding in L 2 (Q, a). 

1.1. Definition. A cubature formula for a finite dimensional space T of functions in 
L 2 (Q, cr) is a pair (X, W), where X is a finite subset of Vt and W : X — ► M*^_ is a weight 
function, such that 

W{x)ip{x) = / if{u)da{uj) 
xex ^ n 

for all <p e T . If W : x \ — ► c(O)/ is the uniform weight, the set X is called a design 
for J 7 . 



x Not functions modulo equality almost everywhere! 
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This definition would make sense for T C L x (0, cr). But the examples we have in mind 
are in L 2 (fi, cr), and we will use Hilbert space methods for which we have to assume that 
T C L 2 (Q, cr); two functions /i, fi G T have a scalar product given by 

(/i I h) = [ fi(oj)f 2 (oj)da(oj). 
Jo. 

There are other denominations for these and related notions, including "Chebyshev qua- 
drature rule" (for our "designs") with appropriate "algebraic degree of exactness"; see, 
e.g., [Gauts-76]. 

We write (• | ■) to denote both the scalar product of two vectors in a Hilbert space of 
functions and the scalar product in the standard Euclidean space W 1 . 

1.2. Example (Simpson's formula). In the case of an interval O = [a, 6] of the real line 

(a < b) and Lebesgue measure, Simpson's formula 

b — a , \ 4(6 — a) ,a + b^ b — a .,. f b , . , 
+ "^6 ^(-2~) + —V{b) = J <f(t)dt 

provides a cubature formula for the space J-^\a,b) of polynomial functions of degree at 
most 3 on [a, b]. More generally, if <p is of class C 4 on [a, 6], we have 
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b — a , x 4(6 — a) /a + 6x, 6 — a /" b . . , 

<fi{a)+ V 1 V (— ) + -g-y(&) ~ / 



< sup 

2880 a <t<6 



Simpson's formula is most often used in its compound form, namely on the n subintervals 
[a + ^-(6 — a) , a + ^(b — a)], 1 < j < n, for n large enough; in this form, it provides a 
very efficient tool in numerical analysis. 



The following notion borrows some of the ingredients of the notion of "polynomial space" 
from [Godsi-93]. 

1.3. Definition. A sequence of polynomial spaces on a finite measure space (O, cr) 
is a nested sequence of finite dimensional spaces 

R = c C • • • C F {k) C • • • 

of functions in L 2 (0, cr) with the following property: J 7 ^ is linearly generated by products 
cpi</?2 with (f \ G and (p2 G jF( fc_1 ), for all k > 1. 



1.4. Remarks, (i) For (F^)k>o as in the definition, there is a natural mapping from the 
fcth symmetric power of onto If n+1 = dim M (J^( 1 )), it follows that dim M (^( fc )) < 

m 

(ii) Many examples of sequences of polynomial spaces have one more property: the 
union 
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1.5. Definition. Let (fi, a) and {F^) k>o be as in Definition 1.3. A cubature formula 

of strength k on O is a pair (X, W), consisting of a finite subset X C O and a weight 
W : X — > which is a cubature formula for J 7 ^. In case W(x) = cr(ft)/ \X\ for all 
x G X, the set X is a geometrical 2 /c-design on fi. 

1.6. Organising questions, (i) For standard examples of finite subsets X of fi, compute 
the largest strength k for which X is a geometrical /c-design. 

(ii) For k given, find designs and cubature formulas (X,W) of strength k with |X| 
minimal. (More on this in Item 1.17, on tight spherical designs.) In case Q is moreover a 
metric space, describe the distance set 

Dx = {c G ]0, oof | c = d(x, y) for some x, y G X, x ^ y} . 

(hi) Dually, for given N, find designs and cubature formulas (X, W) with \X\ = N of 
maximum strength. 

(iv) Asymptotics and equidistribution. Find sequences (X( k \wW) k>1 , where each 
(X^ k \W^) is a cubature formula of strength k, such that the sequence of measures 
( a ^) k>i conver g es to a, where = Ylxexw W^ k \x)8 x and where 8 X denotes the 
Dirac measure of support {x}. Optimize in some sense the speed of convergence. 

When each X^ is a geometrical /c-design, this is the standard question of cx-equidistri- 
bution of a sequence of finite subsets of a measure space; see [Serre-97]. 

In the following proposition, £ 2 (X, W) stands for the Hilbert space of real- valued func- 
tions on X, with scalar product defined by (ipi \ V'2) = Ylxex W (x)^i(x)^2(x) . 

1.7. Proposition. Let (O, a) be a finite measure space. Let H be a finite dimensional 
space of functions in L 2 (0,a). Assume that there exists a finite dimensional subspace TV 
o/L 2 (0, a) such that Ti is linearly generated by products </?iy?2 with cpi, cpz G TV . Let X be 
a finite subset of O and let W : X — > W* + be a weight. 

Then (X, W) is a cubature formula for Ti if and only if the restriction mapping 

p :TV — > £ 2 (X, W) 
is an isometry. In particular, if (X, W) is a cubature formula for Ti, then 
(*) \X\ > dim K (7Y / )- 

Proof. The condition for p to be an isometry is 

^ W(x)(p 1 (x)(p 2 (x) = / (p\{u>)(p2{u>)da{u>) 

for all ipi, (f2 G TV. By hypothesis on Ti and TV, this condition is equivalent to 



xex 

for all (pen. n 



W(x)(p(x) = / ip(oj)da(oj) 



2 In particular cases, we replace "geometrical" by more suggestive adjectives such as "interval" in 
the situation of Example 1.10 with Q an interval in M, "spherical" in the situation of ^ = § n_1 , or 
"Gaussian" if £1 is a real vector space and a a Gaussian measure. However, it should be kept in mind 
that the "Euclidean designs" of [NeuSe-88] and [DelSe-89] (among others) are of a different nature, since 
the reference measure on the Euclidean space Vl = 1™ depends in this case on the weighted set (X, W) 
through some radial factor. 
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1.8. Definition. Let (fi,a) and (J*®) k>o be as in Definition 1.3. A cubature formula 
(X, W) of strength 21 is tight if equality holds in (*) above, with H' = J 7 ^ and H = T^. 

For a spherical cubature formula of strength 21 + 1, "tight" is defined in 1.11. 

Tight cubature formulas are rare: see Discussion 1.17. 

1.9. Quadrature formulas on intervals. Let O be a real interval (a, 6), with — oo < a < 
b < oo (the interval can be open or closed). The spaces J-( k \a, b) of polynomial functions 
on O of degree at most k provide the canonical example for Definition 1.3. 

Consider in particular the case O = [—1,1], with Lebesgue measure, and an integer 
I > 1. Let first 

X® = {x 1 ,i,x 2 ,i,...,xi,i} C [-1,1] 

be an arbitrary subset of / distinct points in [—1, 1]. Let Lj(t) = Yl J.^-x'- } (product over 
% G {1, ...,/}, i 7^ j) be the corresponding elementary Lagrange interpolation polynomials. 
Then (X®, W®) is a "cubature formula" for ^('-^(-l, 1), with the "weight" defined 
on X® by 

W^{x hl ) = J ^L 3 {t)dt (l<j<l). 

But quotation marks are in order since the "weight" values need not be positive; for 
example, in case x^\ = —1 + 2|£j, these Newton-Cotes formulas have all weights 
positive if and only if either / < 8 or / = 10 (see, e.g., Chapter 2 of [CroMi-84] ) . 

Let now X^ be precisely the set of the roots of the Legendre polynomial of degree /, 
namely of P t (t) = (-l)^^r ((l-t 2 ) 1 ). The pair (X^ l \\V^) is a cubature formula 
for JF^-^-M), known as a Gauss quadrature 3 (or sometimes Gauss-Jacobi mecha- 
nical quadrature); this is the unique cubature formula with < / points for JF( 2 ' -1 )(— 1, 1) 
and Lebesgue measure. The weights W^ l \xj : i) of the Gauss formula are strictly positive; 
indeed, since the polynomial Lj defined above have now their squares in JF( 2 ' _1 )(— 1, 1), 
we have 

W®(*i,i) = ]TV^(£ M )L,(£M) 2 = / Lj(t) 2 dt > 
k=i J ~ x 

for j e {1, . . . , /}. (See, e.g., [Gauts-97].) 

In general, there are analogous formulas for other intervals of the real line and other 
measures with finite moments and with dim]R (L 2 ((a, 6), a)) = oo. This is a part of the 
theory of orthogonal polynomials; see Section 3.4 in [Szego-39]. There are also related 
results in larger dimensions; see Section 3.7 in [DunXu-01]. 

1.10. Interval designs. It is a particular case of a theorem of Seymour and Zaslavsky that 
interval designs (named "averaging sets" by these authors) exist for any finite-dimensional 
space of continuous functions on O = ]0, 1[, with various measures; see [SeyZa-84] and 
[Arias-88]. 

An interval design for the space J 7 ^ 21-1 '^— 1,1) of polynomial functions of degree at 
most 21 — 1, and Lebesgue measure, requires strictly more than |/ 2 points (this is a 1937 
result of Bernstein for which we refer to [Krylo-62] and [Korev-94]), and there exists a 

3 If f2 is an interval of M, the term "quadrature" is often used instead of "cubature" as in Definition 1.1. 
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constant c for which it is known that ck 3 points suffice (Section 3.3 in [Korev-94], improving 
on Theorem 2 in [Bajn-91a]). For low values of k, as an answer to a question of Chebyshev 
(1874), Bernstein has shown that there exist interval designs for jF( fc )(— 1, 1) with \X\ = k 
if and only if k < 7 or k = 9, and there is an explicit construction for these values of k (see 
Section 10.3 in [Krylo-62]). 

1.11. Cubature formulas on spheres and tightness. Let O = § n_1 be the unit sphere 
in the Euclidean space M. n , n > 2, and let a denote the probability measure on § n_1 which 
is invariant by the orthogonal group O(n). 

For spheres, the standard example for Definition 1.3 is given by 

R = jr(o)(gn-i^ c jr(i)(§n-i) c ••• C JF^S" -1 ) c ••• 

where JF^^S"- -1 ) denotes the space of restrictions to the sphere of polynomial functions 
on R n of degree at most k. In this paper, when Definition 1.8 is particularized to spherical 
cubature formulas and designs, it is always with respect to this sequence of polynomial 
spaces; thus, a cubature formula (X, W) or strength 21 is tight if 

Consider the space p( fc )(8 n-1 ) of restrictions to S n_1 of polynomial functions on IR n 
which are homogeneous of degree k, and the space 

w (fc)(§»-i) f res trictions to S 71 " 1 of 
polynomial functions on IR n which are homogeneous of degree k and harmonic 4 . We have 

[fc/2] 

p(fc)(gn-l) = w (fc-2j)(gn-l) 
3=0 

for all k > (see e.g. Section IV. 2 in [SteWe-71]). For reasons related to Proposition 3.2 
and Remark 3.3 below, a cubature formula (X, W) or strength 21 + 1 is tight if 

\X\ = 2dim M (V {1) (S 71 ' 1 )) = 2 

Another example for Definition 1.3 is given by 

R = p(°)(§"- 1 ) C P (2) (§ n_1 ) C ••• C ^^(S"- 1 ) C ••• 

(observe that P^^S 71 " 1 ) c V {k+2) (§ n_1 ) for all k > 0, since the restriction of xf + - ■ - + xl 
to § n_1 is the constant 1). 

A subset X C S n_1 is called antipodal if —X = X; such a X can always be written 
(non-uniquely) as a disjoint union of some Y inside X and of —Y. A cubature formula 
(X, W) on S 71 " 1 is antipodal if X is antipodal and if W{-x) = W(x) for all x G X. If 
(X, W) is antipodal, Ylxex W(x)ip(x) = J s „_i (p(u)duj — for any homogenous polynomial 



4 A smooth function (p on 1™ is harmonic if X^?=i 7T~T = 0- 



n + l-1 
n-1 
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function ip of odd degree; it follows that (X, W) is a cubature formula for JF^ 1 )^- 1 ) if 
and only if it is a cubature formula for ^^(S" -1 ). 

1.12. On strength values for spherical designs and cubature formulas. The 

following existence result goes back essentially to the solution of Waring's problem by 
Hurwitz and Hilbert (see [Hurwi-08], [Hilbe-09], page 722 in [Dicks-19], [Ellis-71], and 
[Natha-96]): 

for any I > 0, there exists a cubature formula [Y^ l \ W^) for 7 ? ^ 2 ^(S n_1 ) with 

(^l" 1 ) =dim K (p«(§"- 1 )) < |y« < 

dim M (^ 2 ') = ( n ^ 2/ ~ 1 )-l- 

The lower bound is that of Proposition 1.7; the upper bound is a particular case of Theo- 
rem 2.8 below, and improves by 1 the bound of Theorem 2.8 in [LyuVa-93]. The relation 
with Waring's problem is alluded to in Item 3.9; for the application to Waring's problem, 
it is necessary for the weights of the cubature formula to be rational positive numbers; our 
proof of Theorem 2.8 does not show this, and we refer to Chapter 3 of [Natha-96] for a 
complete proof. 

It follows that there exists an antipodal cubature formula (X"\ W™) of strength 21 + 1 
with 

'n + l-V 



n-1 



< 



<2(" + 2Z - 1 |-2. 
n — 1 



However, the proof of the existence of (Y^ l \ W"') does not provide nice constructions. 
For n > 3 and k > 4, we do not know many good explicit cubature formulas of strength k 



on S n_1 . 



Note the relevant asymptotics. For n fixed: 

'n + l-l\ Z n_1 



n — 1 ) n\ 
For / fixed: 

'n + I — 1\ n l 



n-1 l\ 



if I — > oo. 



if n — > oo. 



The lower bound \X^\ > 2( n ^ 1 1 ) has been improved in some cases, in particular for n 
fixed and / — > oo [Yudin-97]. 

The general result of Seymour and Zaslavsky already quoted in Item 1.10 implies: 
for any n>2 and k > 0, there exists a spherical k-design X C § n_1 . 
The proof of [SeyZa-84] is not constructive and does not give any bound on the size \X\. 

Let X be a non-empty finite subset of S n_1 and let W be a weight on X such that 
E xe xW(x) = l. 

(i) The pair [X, W) is always a cubature formula of strength 0. It is a cubature formula 
of strength 1 if and only if the weighted barycentre J2 x ex W(x)x of (X, W) coincides with 
the origin of IR n . 
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(ii) If (X, W) is a cubature formula of strength 2, then X generates M n . Indeed, if there 
exists a non-zero vector a orthogonal to X, the function u i — ► (a \ u) 2 has a non-zero 
integral on S n_1 but vanishes identically on X. 

Exercise. Denote by Z the Gram matrix of X, defined by Z XjV = (x | y) for x,y G X, 
and by J the matrix with rows and columns indexed by X and with all entries 1. Show 
that X is a spherical 2-design if and only if the three following conditions hold: (ii a ) 
Z XjX = 1 for all x G X, (ii b ) ZJ = 0, (ii c ) Z 2 = n~ l \X\ Z. (Solution in Lemma 13.6.1 of 
[Godsi-93].) 

(iii) If (X, W) is a cubature formula of strength 4, then X cannot be a disjoint union of 
two orthogonal sets; more generally, in case there exist two vectors a, (3 G S n_1 such that 
X C a 1 - U /?-*-, the function cj i — > (a | uj} 2 (/3 | cu) 2 has a non-zero integral on § n_1 but 
vanishes identically on X. 

(iv) If (X, W) is a cubature formula on § n_1 of strength 21 for some I > 1, then the set 
PW(X) = {cj — > (ui | ^) Z } x6X linearly generates "pW(§ n_1 ) . Indeed, assume that some 
p G ^(S™" 1 ) is orthogonal to P®(X). If G ^(S™" 1 ) is defined by 

0(u) = / (w | u) l y(u))dcr{u)) V«el n , 
is™- 1 

the hypothesis means that the restriction of ip to X vanishes; this implies that ip = by 
Proposition 1.7. But </? = implies p = since {c<j — >■ (a; | u ) } u€ §™-i linearly generates 
p(')(S n_1 ). With a terminology borrowed from the theory of lattices (see Reminder 1.15), 
the special case of Property (iv) for / = 2 states that cubature formulas of strength 4 on 
gn-i p rov i(j e perfect sets. 

For cubature formulas on spheres, the following equivalences are useful. Since many 
natural examples are provided by intersections of lattices with spheres of varius radii 
(see 1.16), we find it useful to consider for each p > the sphere p§ n_1 of radius p and 
centre the origin in R n ; we denote again by a the (9(n)-invariant probability measure on 

1.13. Proposition. Consider integers n > 2 and k > 0, a positive number p, a non-empty 
finite subset X of and a weight W : X — > R* + such that XLex W(x) = 1. The 

following conditions are equivalent 5 . 

0) T, x exW(x)<p(x) = f pSn ^p(u)da(u) for all p G V^(R n ) and j G {0,1,..., fc}, 
namely (X, W) is a cubature formula of strength k on pS n_1 . 

(ii) Exex W(x)p(x) = £ xeX W{x)p{gx) for all p G V^(R n ), j G {0, 1, . . . , k}, and 
9 e 0(n). 

(iii) Exex W(x)p(x) = /or all p> G ft^(ffi n ) and j G {1, . . . , fc}. 

(iv) // / zs it/ie largest integer such that 21 < k, there exists a constant C21 such that 

(iv) a w ( x )( x I U ) 2L = P 2l (u I u) 1 V u G M n , 

x6X 



5 Even if there are canonical identifications between V^(R n ), 7>( fe ) (§ n_1 ), and (pS™" 1 ), see Ex- 
ample 1.11, we choose here the first notation. 
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and, if I' is the largest integer such that 21' + 1 < k, then 

(iv) b W(x){x | u) 2l ' +1 = VtiGr. 

xex 

(v) Condition (iv) a holds for any even integer 21 < k and Condition (iv)b holds for 
any odd integer 2V + 1 < k. 

Moreover, the constant in (iv) a is 

1 x 3 x 5 x ••• x {21- 1) 
° 21 ~ n(n + 2)(n + 4) ••• (n + 2/ - 2) 

if I > 0, and c$ = 1 if I = 0. 

Proof. The proof of Theorem 3.2 in [VenMa-01], for spherical designs, readily carries over 
to cubature formulas. □ 

Remarks. If (X, W) is antipodal, then Equality (iv)b is automatically satisfied. 

Besides the conditions of Proposition 1.13, there are many other equivalent ones; see 
[DeGoS-77] and [GoeSe-79], as well as Item 2.14. 



1.14. Examples of spherical designs: orbits of finite groups, inductive cons- 
tructions, distance-regular graphs, and contact points of John's ellipsoids. 

Group orbits. 

It is easy to show that any orbit in § n_1 of any irreducible finite subgroup of 0{n) is a 
2-spherical design; see [Banna-79], but the result can already be found in a 1940 paper by 
R. Brauer and H.S.M. Coxeter (Theorem 3.6.6 in [Marti-03] ) . Moreover, a (2/) -spherical 
design which is antipodal is also a (21 + l)-spherical design. In particular, if (e±, . . . , e n ) is 
the canonical basis of IR n , then the set {±ei, . . . , ±e n } is a spherical 3-design on § n_1 of 
cardinality 2n (hence a tight 3-design); the set {±ei ±e2 ± • • • ±e n } (all choices of signs) is 
another spherical 3-design, of cardinality 2 n , and therefore non tight when n > 3; neither 
of these is a spherical 4-design. 

Orbits of groups generated by reflections provide interesting spherical designs. In par- 
ticular irreducible root systems of type A n , D n , E n are spherical 3-designs (see Chapter V, 
§ 6, Number 2 in [Bourb-68]). Moreover, roots systems of type A 2 , D 4 , E 6 , E 7 provide 
spherical 5-designs and the root system of type Eg provides a spherical 7-design. 

For n > 3, there is a maximal strength k max (n) for spherical designs in § n_1 which are 
orbits of finite subgroups of 0(n). It is moreover conjectured that sup n>3 k max (n) < oo; the 
maximal value of k we know is k = 19 and occurs in dimension n = 4 (see 3.8 below). For 
more on spherical designs which are orbits of finite groups, see [Banna-79] and [HarPa-04]. 

Inductive constructions. 

Spherical designs can be constructed inductively as follows. On S 1 , the vertices of a 
regular iV-gon constitute a /c-design if and only if N > k + 1. For n > 3 and some given 
k, assume that we have a spherical /c-design Y on S n_2 , as well as an interval design 
Z C ] — 1, 1[ for the space jF( fc )(— 1, 1) of polynomial functions of degrees at most k and for 
the measure da(t) = (1 - t 2 ) (n - 3) / 2 dt. Then 



x = u x i-\\zr y ,z) gs 11 



y E Y and z £ Z 
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is a spherical /c-design in S n_1 . Denote by M' n (k) the smallest integer such that, for any 
N > M' n (k), there exists a spherical /c-design of size N on §> n_1 . With the construction 
above, it can be shown that 

/ n(n+l) „\ 

M' n {k) < 0(fc^" 2 ). 

It is conjectured that n( - n 2 +1 - ) — 1 can be reduced to n( - n 2 ~ 1 - ) . For all this, see [Bajn-91b] 
and other papers by Bajnok. 

In case n = 3, it is a result of [KorMe-93] that, for all k > 0, there exists a spherical 
/c-design of size 0(k s ). 

Constructions of spherical designs appear also in [Mimur-90] , [Bajno-92], and [Bajno- 
98]. 

Distance regular graphs. 

Consider a distance-regular graph T, with vertex set V(r), of valency k. Let 9 be an 
eigenvalue of T, 9 ^ k, of some multiplicity to; let p denote the orthogonal projection of 
the space of functions V(F) — > R to the #-eigenspace of the adjacency matrix of V. Then 
p(V(T)) is a spherical 2-design in pS m_1 for the appropriate radius p (Corollary 13.6.2 in 
[Godsi-93]). 

Contact points of John's ellipsoids. 

Let K be a convex body in n-space such that the maximal volume ellipsoid in K is the 
unit n-ball B. It is a theorem of John that there exist a finite subset X of K n dB and a 
positive weight function W on X such that (X, W) is a cubature formula for dB = § n_1 
of strength 2, indeed of strength 3 if K is assumed to be symmetric; see Lecture 3 in [Ball- 
97] . It could be interesting to investigate systematically polyhedra K such that K n dB 
is a spherical /c-design, or such that there exists a weight W for which (K n dB, W) is a 
cubature formula of strength k, for various values of k. 

1.15. A reminder on lattices. Let V be a Euclidean space of dimension n > 1. A 
lattice is a discrete subgroup A of V which generates V as a vector space, and n is the 
rank of A. If A C V is a lattice, so is its dual A* = {x G V \ (x \ A) C Z}. A lattice 
A is integral if (A | A) C Z, equivalently if A C A*. A lattice A is unimodular if 
Vol(V/A) = 1; an integral lattice is unimodular if and only if A = A*. A lattice A is even 
if (A | A) G 2Z for every A G A; an odd lattice is an integral lattice which is not even. 
Two lattices A C V, A' C V are equivalent if there exists an isometry g from V onto V 
such that g(A) = A'. There is a natural notion of orthogonal direct sum of lattices, and a 
lattice is irreducible if it is not equivalent to a direct sum A © A' with A, A' ^ 0. 
Let A C V be an integral lattice; for to > 1, we denote by 

A m = {A G A | (A | A) = to} 

the shell of radius to (namely of radius ^fm in the usual sense of Euclidean geometry). 
For any A G A l5 there is an integral lattice A' in the orthogonal subspace A^ such that 
A = ZA © A'; in particular, Ai = for an irreducible lattice of rank at least 2. 

A standard example of an odd unimodular lattice is the cubical lattice Z n C M. n . If 
n < 11, any odd unimodular lattice is equivalent to one of these. 
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Let n be a multiple of 4. In the standard Euclidean space M n , define the lattice D n = 
{A G Z n | EF=i A * = ( mod 2 )} and the Witt latticed = D„u((|, . . . , J) + £> n ). Then 
r n is integral and unimodular; moreover, it is even if and only if n is a multiple of 8. 
Even unimodular lattices exist only in dimensions n = (mod 8). If n = 8, any even 
unimodular lattice is equivalent to Eg, also known as the root lattice of type E% or 
the Korkine-Zolotareff lattice. If n = 16, any even unimodular lattice is equivalent to 
either E 16 (which is irreducible) or to E 8 © E 8 . If n = 24, there are 24 equivalence classes 
of even unimodular lattices (Niemeier's classification, 1968); the most famous of them is 
the lattice discovered by Leech (1964), which is the only even unimodular lattice A in 
dimension n < 31 such that A2 = 0, and which has a remarkably high density (see, e.g., 
Table 1.5 in Chapter 1 of [ConSl-99]). If n = 32, the number of equivalence classes of even 
unimodular lattices is larger than 8 x 10 7 . 

There is a classification of integral unimodular lattices of small rank. For n < 16, 
Kneser (1957) has shown that the only irreducible integral unimodular lattices are Z, E n 
for n = 8, 12, 16, and three odd lattices in dimensions 14, 15, and 16. In particular, if 
9 < n < 11, any odd unimodular lattice is equivalent to either Z n or E 8 © Z n ~ 8 . Integral 
unimodular lattices have been later classified for n < 23 (Conway and Sloane, 1982) and 
n < 25 (Borcherds, 1984). For n < 24, see Chapters 16 and 17 of [ConSl-99]; minor 
corrections to previous tables are given in [Bache-97]. 

Let A C V be a lattice (possibly neither integral nor unimodular). For r > 0, let B r 
denote the ball or radius r centred at the origin of V . The sphere packing associated to 
A is Ua6a( A + -^ r )' w here r denotes the largest real number such that the balls A + B r 
have disjoint interiors (A G A); the density of A is the number 



The lattice A is extreme if this density is a local maximum in the space of all lattices of 
dimension n; since density is homothety-invariant, a unimodular lattice is extreme if its 
density is a local maximum in the space of all unimodular lattices of dimension n, a space 
which can be identified with the double coset space SL(n, Z) \ SL(n 7 M)/0(n). Denote by 

Ashort = {/U G A , fj, ^ \ (jjl \ /jl) < {X \ X) for all A e A , A ^ 0} 

the shell of short vectors in A; the lattice A is eutactic if there exists a weight W : 
A s hort — ► such that (A s hort 5 W) is a spherical cubature formula of strength 3, and 
perfect if the set {lv i — > (uj | A) 2 | A G A s hort} linearly generates the space of homogeneous 
polynomial of degree 2 on V. The first result involving both lattices and cubature formulas 
is the following theorem of Voronoi (1908): 



For A to be extreme, it is sufficient that A is strongly perfect, which means that A s h or t 
is a spherical 5-design (Theorem 6.4 in [VenMa-01]). Strongly perfect lattices have been 
classified in dimensions n < 11, where they occur in dimensions n = 1, 2, 4, 6, 7, 8, 10 only 
[VenMa-01], [NebVe-00]; other examples occur in [BacVe-01]. 

Standard references on lattices include [ConSl-99], [Ebeli-94], [Marti-03] , [MilHu-73] , 
[Serre-70], and [VenMa-01]. 

1.16. Examples of spherical designs: lattice designs. Let A be a lattice in M n . Any 
non-empty shell A m provides a spherical design of some strength in the sphere y^S" -1 . 



1^ Vol (( |J (A + B r )) n / Vol(S fl ) 



Vol(£ r ) 
Vol(F/A) ' 



a lattice is extreme if and only if it is both eutactic and perfect. 
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This connection goes back to [Venko-84]. See [Pache] for many examples; let us indicate 
here a sample. 

If A C M 8 is a root lattice of type Eg, then is a spherical 7-design for any m > 1. 
Moreover, it can be shown that a shell A2 m is a 8-design if and only if the Ramanujan 
coefficient r(m) of the modular form 

oo oo 

A 24 (*) = q 2 X\{l-q 2m ? A = $>(m)<z 2m 

m=l m=l 

= q 2 — 24q 4 + 256q 6 - U72q 8 + ■■■ (q = e mz ) 

is zero. (See Example 3.10 below for a second appearance of this modular form A24 of 
weight 12.) Now it is a famous conjecture of Lehmer [Lehme-47] that r(m) 7^ for all 
m > 1. If T C 1R 16 is a Witt lattice (an irreducible even unimodular lattice, uniquely 
defined up to isometry by these properties in dimension 16), then ]?2 m is a spherical 3- 
design for any m > 1, and the condition for one of these shells to be a 4-design happens 
to be again the vanishing of the corresponding Ramanujan coefficients. The same holds 
for the reducible even unimodular lattice A® A 6 IR 16 . Consequently, the following four 
claims are simultaneously true or not true: 

(i) Lehmer 7 s conjecture holds, namely r(m) 7^ for all m > 1; 

(ii) no shell of the root lattice A of type E$ is a spherical 8-design; 
(hi) no shell of the lattice A © A e ffi 16 is a spherical 4-design; 

(iv) no shell of the Witt lattice V e R 16 is a spherical 4-design. 

Claim (i) has been checked for m < 10 15 [Serre—85]. 

It is easy to formulate other equivalences of the same kind. For example, any shell of the 
Leech lattice L (see 1.15) is a spherical 11-design, and the condition for £2™, (m > 2) to 
be a 12-design is that the mth coefficient ju(m) = Yl^i T (j) T ( m ~ j) °f (A24) 2 vanishes. 
Thus 

A*m 7^ for any m > 2 if and only if Li2m is not a spherical 12-design 
and, moreover, we have checked that \i m ^ for m < 1200. We do not know of any other 
lattice with a shell which is a spherical design of strength > 11. 

Similarly, let now A denote a Niemeier lattice, namely an even unimodular lattice in 
dimension 24 with A 2 7^ (up to isometry, there are 23 such lattices). Denote by Q the 
modular form of weight 4 defined by 

Q(q) = (Mi)f - -^(Mq)h(q))' 

= (E^) 8 -( E ^y(£(-«r 2 ) 4 

= 1 + 240q 2 + 2160/ + • • • 

and set 

00 

Q(q)A 24 (q) = »mq m . 

m=l 

Then 

v m 7^ if and only if A 2m is not a spherical 4-design (m > 1) 
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and, moreover, we have checked that v m ^ for m < 1200. (Observe that the condition 
is the same for any of the 23 Niemeier lattices.) 

For one more example in this class, consider the cubical lattice Z n . The criterion stated 
in 1.14 in terms of irreducible representations of finite subgroups of 0(n) shows that all 
non-empty shells are spherical 3-designs. Moreover, it can be shown that there are two 
classes of "special shells" which are spherical 5-designs, namely (Z 4 ) m for m = 2a and 
(Z 7 ) m for m = 4 b (8a + 3). Let us restrict now for brevity to n > 8, and denote by 
the modular form of weight 4 + n/2 defined by 

1 oo 

® [n] (q) = -(e 2 (q)eMf (OMY = E* n - 

m=l 

Then 

Km 7^ if and only if ( r L n ) m is not a spherical 4-design (m > 1) 
and, moreover, we have checked that these hold for all n > 8 and m < 1200. (See [Pache], 
which contains a discussion including 1 < n < 7.) 

1.17. Tight spherical designs on spheres of dimension n — 1 > 2. (See [BaMuV-04] 
and references there.) 

Even strengths. 

Tight spherical (2/)-designs do not exist when 21 > 6. 

Tight spherical 4-designs in S™" 1 are of size dimR^ 2 ^" -1 )) = \n(n + 3). They 
cannot exist unless n is of the form (2m + 1) 2 — 3. If m = 1 or m = 2 examples are known, 
and known to be unique up to isometry; they are respectively of size 27 in § 5 and size 275 
in § 21 (see below). If m = 3 and m = 4 (and infinitely many larger values), non-existence 
has been proved. 

Tight spherical 2-designs exist in all dimensions, and are regular simplices. 
Odd strengths. 

For I > 0, a tight spherical (21 + l)-design is necessarily antipodal, by Theorem 5.12 in 
[DeGoS-77]. 

Tight spherical (21 + l)-designs do not exist when 2/ + 1 > 9, up to one exception (which 
is unique up to isometry): the 196 560 short vectors of a Leech lattice which provide (after 
dividing all vectors by 2) an 11-design in § 23 . Observe that 



196560 = 




= 2dim M (p( 5 )(§ 23 )). 



Tight spherical 7-designs in § n_1 are of size |n(n + l)(n + 2). They cannot exist unless 
n is of the form 3m 2 — 4. If m = 2 or m = 3, examples are known, and known to be 
unique up to isometry. (Up to homothety, they are respectively a root system of type E$ 
and the short vertices of the unimodular integral lattice denoted by O23 in [VenMa-01] 6 .) 
If m = 4 and m = 5 (and infinitely many larger values), non-existence has been proved. 

6 Let L be a Leech lattice. Let first e 6 L be a short vector, with (e | e) = 4. Denote by p the 
orthogonal projection of R 24 onto e 1 - and set L' e = {x £ L \ (e \ x) = (mod 2)}. Then O23 = p(L' e ) is a 
unimodular integral lattice with min{(a; | x) \ x G O23, x 7^ 0} = 3 and \{x 6 O23 | (x \ x) = 3}| = 4600. 
Let then / 6 L be a vector such that (/ | /) = 6. Then M23 = Lfl f 1 - is an integral lattice with 
min{(a; | x) \ x € M23, x 7^ 0} = 4; if M| 3 denote its dual lattice, then [M| 3 : M23] = 4, and there are 552 
short vectors in MX„. 
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Tight spherical 5-designs in § n_1 are of size n(n + 1). They cannot exist unless n is 
either 3 or of the form (2m + l) 2 — 2. If n = 3, or m = 1, or m = 2, examples are known, 
and known to be unique up to isometry. (They are respectively a regular icosahedron and, 
up to homothety, the short vectors of a lattice which is dual to a root lattice of type E7 
and the short vectors of a lattice of type M| 3 , with the notation of [VenMa-01].) If m = 3 
and m = 4 (and infinitely many larger values), non-existence has been proved. 

Let X C be a tight spherical 5-design. It is a consequence of Proposition 3.4 

that (x I y) = ±a for any x,y G X with x 7^ ±y, where a = 1/v^n + 2. Choose e E X 
and set X$ = {x G X | (x | e) = a}. Then Xo is a tight spherical 4-design in the sphere 
Vl — a 2 S n_2 centred at ae in the affine hyperplane {x G M n | (x | e) = a}, by Theorem 8.2 
in [DeGoS-77]. 

Tight spherical 3-designs exist in all dimensions, and are of the form {±ei, . . . , ±e n }, 
where {ei, . . . , e n } is an orthonormal basis of M n . Tight spherical 1-designs are of the form 
{±ei}- 

For results on tight designs on other spaces (0,cr), see [Babal] and [Baba2]. 

1.18. Designs with few points on § 2 . 

Given an integer N > 1, it is a natural question to ask what is the largest integer fcjv > 
for which there exists a spherical /cTv-design with iV points on a sphere, say here on S 2 . 
We report now some answers to this question, from [HarSl-96]. 



N = 1 


h = o 


single point 


N = 2 


k 2 = 1 


pair of antipodal points 


N = 3 


fe 3 = l 


equatorial equilateral triangle 


N = 4 


k 4 = 2 


regular tetrahedron (tight) 


N = 5 


&5 = 1 




N = 6 


k 6 = 3 


regular octahedron (tight) 


N = 7 


fc 7 = 2 




N = 8 


fc 8 = 3 


cube (non-tight) 


N = 9 


kg = 2 




N = 10 


feio = 3 




N = 11 


hi = 3 




iV = 12 


fel2 = 5 


regular icosahedron (tight) 


< iV < 23 


A:tv e {3,4,5} 




iV = 20 


&20 = 5 


regular dodecahedron (non tight) 


iV = 24 


/c 2 4 = 7 


improved snub cube 


N = 60 


/C60 = 10 





Remarks, (i) The table shows that /cat is not monotonic as a function of N. 

(ii) We know from Propositions 1.7 and 3.2 that, on S 2 , any spherical (2/)-design has 
size N > (I + l) 2 and any spherical (21 + l)-design has size N > (I + 2){l + 1). The table 
above shows that these bounds are not sharp unless 21 G {0, 2} or 21 + 1 G {1, 3, 5}. For 
example, a 4-design has minimal size 12(> 9) and a 7-design has minimal size 24(> 20). 

(iii) The table refers to isolated examples; here is a continuous family, from [HarSl-92]. 
Distribute the 12 points of a regular icosahedron into two poles iV and S, and two sets 
P, Q of 5 points each in planes orthogonal to the diameter joining N to S. Let Xq denote 
the union of the poles, the set P, and the image of the set Q by a rotation of angle 9 fixing 
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N and S. If 9 is not a multiple of 27r/5, then Xq is a spherical 4-design which is not a 
5-design. 

(iv) The Archimedes' regular snub cube 7 (24 vertices) is a spherical 3-design (not a 
4-design), while the improved snub cube reported to in the table is indeed a spherical 
7-design. 

(v) The regular truncated icosahedron (= soccer ball, 60 vertices) is a spherical 5- 
design (not a 6-design). The "improved soccer ball" of [GoSe-81a/b], which is almost 
indistinguishable from the regular one with the naked eye, is a 9-design. The spherical 
10-design of size 60 which appears in [HarSl-96] is quite different. 

(vi) The same paper shows a 9-design with N = 48 points; indeed k^ = 9 and k^ < 9 
whenever iV < 48 or iV G {49, 51, 53}. On the other hand, there exists a cubature formula 
on § 2 of strength 9 and size 32 (Section 5 in [GoSe-81a], and Item 3.7 below). 

(vii) It is conceivable that many values of iV are relevant for extra-mathematical reasons: 
N < 12 pores on pollen-grains (botany); 

N > 60 atoms in various large carbon molecules (chemistry); 

N ~ 20 000 detectors for a PET tomography of the brain (medical imaging) ; 

100 < N < 10 20 charged particles on a conducting sphere (electrostatics). 

In Section 3 below, there are other examples of cubature formulas and designs on spheres. 

1.19. Several quality criteria for spherical configurations. What is the best way to 
arrange a given number N of points on a sphere, say here on S 2 ? The answer depends of 
course of what is meant by "the best" . Besides maximizing the strength of the configuration 
viewed as a spherical design, we mention here two other criteria. 

The Tammes' problem asks what are the configurations {x\, . . . ,x n } C S 2 which 
maximize m.mi<i j< n (distance^, Xj)), or equivalently minimize maxi< i j< ni j^j((a;i | 
Xj)). Configurations of N points on § 2 are discussed in [Tamme-28] for N < 12. Pieter 
L.M. Tammes (1903-1980) is a Dutch botanist who was interested in the distribution of 
places of exit on pollen grains (the most frequent case seems to be N = 3); he should not 
be confused with his aunt Tine Tammes (1871-1947), who made important contributions 
to early genetics [Stamh-95]. 

Here are some of the configurations which are the best from the point of view of Tammes' 
problem: those of Example 1.18 for N = 4, 6, 12 (regular polytopes with triangular faces), 
but other configurations for N = 8 (square antiprisms), N = 20 (unknown configurations 
which are not regular dodecahedras [vdWae-52]), and iV = 24 (snub cubes, see ([SchWa- 
51] and [Robin-61]). More on Tammes' problem in mathematics in [Coxet-62], in Section 
35 of [Fejes-64], in Section 2.3 of Chapter 1 of [ConSl-99], and in Chapter 3 of [EriZi-01]. 

Let X be a finite subset of § n_1 . When the emphasis is on properties like those of 
Proposition 1.13, the set X is called a spherical design. However, when the emphasis is on 
the distance set Dx defined in 1.6, the standard name for X is that of spherical code. 
For many constructions of spherical codes, see [EriZi-01]. 

Problem 7 of Smale's Mathematical problems for the next century [Smale-98] asks what 
are the configurations which maximize rii<i<j<jv distance^, Xj). The problem is mo- 
tivated by complexity theory, and the search of algorithms related to the fundamental 
theorem of algebra. 



See for example 



http://mathworld.wolfram.com/SnubCube.html 
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For a discussion of related criteria, see [SaaKu-97] . 

1.20. Cubature formulas on other spaces. There are documented cubature formulas 
of strength 3,5, .. . (with respect to the space of polynomials) on the hypercube [—1, l] n 
of M 71 . For compact subsets of the plane with positive area, there are cubature formulas of 
strength k and of size \{k + l)(k + 2). See [DavRa-84], in particular Section 5.7. 

Let O = T 2 be a 2-torus of revolution embedded in the Euclidean space 1R 3 , together 
with its standard probability measure a (up to a normalization factor, a is the area given 
by the first fundamental form of the surface T 2 embedded in 3-space). Let jF( fc )(T 2 ) be the 
space of functions on the torus which are restrictions of polynomial functions of total degree 
at most k on M 3 . Kuijlaars [Kuijl-95] has shown that there exist constants C\,C2 > 
with the following property: for any k > 0, there exists a geometrical design on .^(^(T 2 ) 
with a number N of points satisfying C\k 2 < N < C 2 k 2 . 

2. Cubature formulas and reproducing kernels 

This section begins with a reminder of standard material which goes back at least to 
[Arons-50] and [Krein-63]; see also [Stewa-76] and [BeHaG-03]. 

Let 7i be a real Hilbert space of functions on a set Q. We denote by (ipi | ^2) the scalar 
product of if 1 and if 2 in Ti- We assume 8 that, for each u G O, the evaluation H — ► R, 
tp 1 — ► <p(w) is continuous; this implies that there exists a function ip u G T~L such that 
(p(uj) = (if I (fij) for all (p G TL. 

2.1. Definition. With the notation above, the reproducing kernel ofTC is the function 
1:0x0 — > K. defined by ^(u',u) = (</? w | <^v) for all G O, or equivalently by 

The terminology, "reproducing", refers to the equality <p(u) = (<£>(•) | ^ , (-,^)) for all 
(f G H and w6 0. 

2.2. Proposition. Let fi, (<A»>) w gn> an d ^ ^ e a5 above. 

(i) The kernel $ is of positive type. In particular, its diagonal values $(u;, a;) are positive, 
and |$(u;, u/)| 2 < $(cu, u;)<I>(u/, a;') /or a// a;, a;' G O. 
fii,) The family (v 7 w) a , 6 n generates 7i. 
(Hi) IfTC 7^ 0, t/ien <&(a;,u;) 7^ /or some wGO. 
firj If {ej)- eJ is any orthonormal basis ofTL, then 

jeJ 

for all u>, u' G O. 

(i>,) In case there exists a finite positive measure a on Q, such that TL is a closed subspace 
o/L 2 (fi,a), with J nxCl \^(u,u')\ 2 da(u)da(u') < 00, then 

diniR(7Y) = / &(uj,uj)da(uj) < 00. 
Jn 



8 In all examples appearing below, TL is finite dimensional and this condition is therefore automatically 
fulfilled. Hence, the reader can assume to start with that dim.^(TL) < 00. 
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(vi) Assume moreover that Q is a topological space and that the dimension ofTL is finite. 
Then the following conditions are equivalent: 
(vi a ) all functions in 7i are continuous; 
(vi ) the kernel $ : Q x Q — > R is continuous; 
(vi c ) the mapping uj i — > (p u from O to 7i is continuous. 

Proof, (i) By definition, the kernel $ is of positive type if it is symmetric and if 

n 

(*) Yl X 3 X k$(Uj,UJ k ) > 

j,k=l 

for all integers n, real numbers Ai, . . . , A n and points u>i, . . . , oj n in O. This is clear here, 
since the left-hand term Y^Jj k=i x j x k (<£uj k \ <fu>j ) of (*) is equal to the square of the Hilbert- 
space norm of the sum Ylk=i x kfuj k - The last claim of (i) follows by the Cauchy-Schwarz 
inequality. 

(ii) Observe that, for <p G TL, the condition (<p \ Lp u ) = for all u G O implies that 
<fi(u>) = for all w G O. 

(Hi) Assume that $(oj, uj) = for all wfSl. Then, by (z), <E>(u;, uj') = for all uj, uj' G O, 
and it follows from (ii) that H = 0. 

(iv) Evaluate the Fourier expansion (p u i = Ej € j(¥V I e i) e i a ^ ^ ne P°i n t ^ 5 an d obtain 
$(w,a/) = Ei6j(<A.' I ej)(ej | p w ) = E je j e jM e jV)- 

(v) Denote by K$, : L 2 (0,cr) — >■ L 2 (0,cr) the linear operator defined by the kernel $, 
namely by (K$(p)(u) = j Q &(u, ui')(p(ui')da(uj') for all <p G Ti. and w£fl. On the one hand, 
K$> is a Hilbert- Schmidt operator, because of the L 2 -condition on on the other hand, 
K$> is the identity, since the kernel $ is reproducing. It follows that the dimension of Ti is 
finite. Moreover, we have 

/ Q(uj,uj)do~(uj) = / \ej(uj)\ 2 da(uj) = ^^||ej|| 2 = dim^(H) 
Jq jeJ^ Q jeJ 

by (iv). 

(vi) Since dim^(T-[) < oo, the mapping of (vi c ) is continuous if and only if the real- 
valued functions u i — ► (<p \ <p u ) = <*p(ui) are continuous for all <p G TC, so that (w ) and 
(ui c ) are equivalent. 

If (vi c ) holds, then $ is continuous since it is the composition of the continuous mapping 
(u>',uj) i — ► (<fu!,<fiu!') from O x O to 7i x with the scalar product. If (vi ) holds, then 
ifuj = $(',oj) is continuous for all w£ll, and (ui a ) follows by (zz). □ 

2.3. Example: atomic measure. Let X be a set and W : X — > a positive-valued 
function. The Hilbert space £ 2 (X, W), with scalar product given by 

(ipi Ifa) = W(x)Mx)Mx) for all Vi, ^2 G £ 2 (X, W), 
xex 

gives rise to the functions 

if) - 1 S ■ y i ( if y = X 
1 ) otherwise 
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and to the reproducing kernel \1/ with values 

y) = (ip y | ip x ) = 



W(x)- 1 if y = x 
otherwise. 



In particular, if X is finite and W(x) = \X\ 1 for all x G X, then \X\ 1 ^ is the charac- 
teristic function of the diagonal in X x X. 

2.4. Example: standard reproducing kernels on spheres. Let the notation be as in 
Example 1.11 and let k > 0. Each of the spaces ^(^(S" -1 ), ^(S™" 1 ), ^(S™- 1 ) is a 
subspace of the Hilbert space L 2 (S> n_1 , a). 

There exists a unique polynomial Q^ k \T) G R[T] of degree fe with the following pro- 
perties: for any u> G S n_1 , the polynomial function defined on M n which is homogeneous 
of degree k and of which the restriction to S n_1 is given by 9 

^ . j g(*)(( w | w ')) 

is harmonic, and Q( fe )(l) = dim ]R (7i( A: )(S n ~ 1 )); the polynomial is a form of a Gegen- 
bauer polynomial (see, e.g., Theorems IV.2.12 and IV. 2. 14 in [SteWe-71], or Section IX. 3 
in [Vilen-68]). It is routine to check that Q W (X) = 1, Q (1) (X) = nX, and 

k + 1 Q( k+1 \X) = XQ^{X) - n + k ~ 3 A Q^ k - 1 \X) 



n + 2k v ' ' n + 2k-4 

for k > 2. Observe that Q^ k \T) is of the form 

[fc/2] 

g( fc )(T) = ^2(-iy c f ) T k - 2 ^ 

j=0 

with Cq \ c[ h \ . . . > 0. The reproducing kernel of the space is given by 

& k \u,u') = Q {k) {{u | u')) 
for all G S n_1 . It is a restriction of the kernel 

[fc/2] 

l n xt n 9 (w,u;') i— > ^(-l) J cJ fc) (a; | u') k ~ 2j ' (u \ u) j (J \ G R 

j=0 

which is homogeneous of degree k in each variable uj,oj' separately. 

The reproducing kernels of the spaces p( fc )(§ n_1 ) and jF( fc )(S n_1 ) are given similarly in 
terms of the polynomials 

[fc/2] k 

C( fc )(T) = Y,Q (k ~ 2j) (T) and R^ k \T) = ^Q^(T). 
i=o j=o 



3 Where (u \ u') = £?=i w^. 
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Observe that 



dim R (^W(S n - 1 )) = QW(1) = 
dim M (p( fc )(§ n - 1 )) = C^(l) = 
dim R (.F (fc) (S n - 1 )) = R( k \l) = 




+ 



( 



( 



n + k-2 
n-1 



n + k — 3 
n — 1 



) 



by Proposition 2. 2. v. 

This material is standard: see, e.g., [SteWe-71] and [DeGoS-77]. 

Proposition 2.6 is the first general existence result for cubature formulas of this exposi- 
tion. (It is a minor strengthening of Theorems 2.8 and 3.17 in [LyuVa-93].) 

2.5. Lemma. Let V be a finite dimensional real vector space, let n be a probability measure 
on V , and let b = f v vdfx(v) denote its barycentre. If C is a convex subset of V such that 
H(C) = 1, then beC. 

Remarks, (i) The lemma is probably well-known, but we haven't been able to trace it in 
print. The proof below was shown to us by Yves Benoist. 

(ii) The point of the lemma is that b is not only in the closure of C, but in C itself. 

Proof. Let U be the minimal affine subspace of V such that fi(CnU) = 1; upon restricting 
C, we can assume that C C U. Upon translating \x and C, we can assume that 6 = 0. 

We claim that b = is in the interior of C inside U. If this was not true, there would 
exist by the Hahn-Banach theorem a non-zero linear form £ on U such that £(c) > for 
all c G C. The equality 



would imply £(c) = for almost every c G C; thus, we would have n(C fl Ker(£)) = 1, in 
contradiction with the definition of U. 

Hence b is in the interior of C inside U, and in particular b G C. □ 

2.6. Proposition. Let (O, a) be a finite measure space. Let H be a finite dimensional 
Hilbert space of functions on O which is a subspace of L 2 (0, a) and which contains the 
constant functions. 

Then there exists a cubature formula (X, W) for H such that 



Proof. Assume first that a is a probability measure. Let H° be the orthogonal supplement 
of the constants in H and let $° denote its reproducing kernel. Recall that, for wefl, the 
function v?° : u' — ► is in H°. The set 




X\ < dim R (ft). 



O = G H° | (p = (p^ for some uj G 0} 
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linearly generates H° by Proposition 2.2.ii. Observe that 

/ ipl(u>')dcr(u>) = [ ^,{uj)da{uj) = (1 | ^,) n0 = 

for all u>' G fi, by definition of H°. It follows from the previous lemma, applied to the 
image of a on fi, that is in the convex hull of O. 

By Caratheodory's theorem (see e.g. Theorem 11.1.8.6 in [Berge-78]), there exist a 
finite subset X of O of cardinality at most d\m R {7i G ) + 1 = dim^iTi) and a weight function 
W : X — ► R* + such that 

W ( x ) = 1 and W ( x ) l fx = G H°. 

xEX xEX 

For all x G X and u' G O, we have <p° x (u') = $°(o/, x) = a;') = </?°/(x); hence 

xex J n 
for all oj' G O. Since ( ( / : '°') w / ef 2 g enera tes H°, we have 

W(x)<p(x) = / <^(o;)cicr(a;) 

for all (/? G H°. Since the same equality holds obviously for constant functions, it holds 
also for all (p G TL 

In case a is not a probability measure, we can multiply all values of W by cr(fi). Then 
the previous equality holds again for constant functions on O, and a posteriori for all 

<p g n. □ 

In the situation of Proposition 2.6, suppose moreover that we have a cubaure formula 
(X', W') for 7Y. Then there exists a cubature formula (X, W) for 7i such that 

X CX' and |X| < dim M (7Y). 

(This follows from the proof in [Berge-78] of Caratheodory's theorem. Alternatively, we 
can apply Proposition 2.6 to the subspace of £ 2 (X\ W') of restrictions to X' of functions 
in TL.) 

2.7. Proposition. In the situation of the previous proposition, assume moreover that Q 
is a connected topological space and that H is a space of continuous functions. 
Then the bound on the size of X can be improved to 

\X\ < dim R (H) - 1. 

Proof. With the notation of the proof of Proposition 2.6, the subspace Q of 7i° is a continu- 
ous image of Q by Proposition 2.2.vi, so that O is connected. In this case, the Caratheodory 
bound on X can be lowered by 1, by a classical theorem of Fenchel. See, e.g., Theorem 18 
in [Eggle-58]. □ 



The following result is an immediate consequence of Propositions 1.7 and 2.7. 
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2.8. Theorem (existence of cubature formulas, with bounds on sizes). Let Q be 

a connected topological space with a finite measure a and let (jF( fc )) fc>Q be a sequence of 

polynomial spaces on (fi, a). Assume that functions in U/eXD^ 7 ^ are continuous on O. 
Choose I > 0. 

Then there exists a finite subset X of O such that 

dim K (jFW) < \ X \ < dimM^ 20 ) - 1. 

and a weight W : X — > R!J_ such that (X, W) is a cubature formula of strength 21 on O. 

Our next proposition (2.10) shows properties of cubature formulas (X, W) for which the 
previous lower bound is an equality, namely which are tight (Definition 1.8). 

2.9. Lemma. Consider data consisting of 

(a) a finite measure space (Q, a), a Hilbert space of functions H C L 2 (fi, a), and the 
corresponding reproducing kernel <£>; 

(b) a finite subset X of Vt, a weight W : X — > the Hilbert space £ 2 (X, W), and 
the kernel ^> as in Example 2.3; 

(c) the restriction mapping p : 7i — > £ 2 (X, W), tp \ — > <f \x, o,nd the adjoint mapping 
p* : £ 2 (X, W) — >n. 

With the previous notation, namely with ip^ as in 2. 1 and ifj x as in 2. 3, we have 

(i) p*(ip x ) = <p x for all x E X; 

(ii) p* is an isometry if and only if ^ is the restriction o/$ to X x X; 

(iii) in case dim^(H) = \X\, the mapping p is an isometry if and only if ^ is the 
restriction o/$ to X x X. 

Proof, (i) By definition of p* , we have for all (p G H and for all x E X 

(ip | p*{ip x )) = (p(<p) I i>x) = (p(<p))(x) = <p(x) 

and therefore p*(t(j x ) = ip x . 

(ii) Let x,y E X. On the one hand, (p*(ip x ) | p*{ip y )) = $>(x,y) by (i); on the other 
hand, (ip x \ ip y ) = ^(x,y). Claim (ii) follows. 

(iii) A linear mapping between two finite-dimensional Hilbert spaces of the same dimen- 
sion is an isometry if and only if its adjoint is an isometry ; thus (iii) follows from (ii). 
□ 

2.10. Proposition. Let (O, a) be a finite measure space and let (J I '^)k>o be a sequence 
of polynomial spaces on (O, a); assume that the corresponding reproducing kernels <£>( fc ) are 
constant on the diagonal ofVt. 

Choose I > and let (X, W) be a cubature formula of strength 21 on Q which is tight. 
Then W is uniform, namely X is a tight geometrical {2l)-design. 

Proof. The restriction mapping 



p:T {V) — >£ 2 (X,W) 
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is an isometry by Proposition 1.7. If (X, W) is tight, the domain and the range of p have 
the same dimension, so that p is onto. Thus, if ^> : X x X — > R denotes as in Example 2.3 
the reproducing kernel of £ 2 (X 7 W), then 

^(x,x') = x') for all x,x'eX 

by Lemma 2.9. It follows that ^ is constant on the diagonal of X x X, so that the weight 
W is constant by Example 2.3; otherwise said, X is a geometrical (2/)-design. □ 

2.11. Definition. Let O be a metric space, let a be a finite measure on Vt, let TL be a 
finite dimensional Hilbert space of functions on O which is a subspace of L 2 (0, a), and let 
$ denote the corresponding reproducing kernel. We say that Condition (M) holds for 
7i and a function S if the values of the kernel depend only on the distances: 

$(u,u}') = S(d(oj,u')) 

for all u>, ui' G Vt, and for some function S : R+ — > R. 
Observe that, by Proposition 2.2.v, this implies that 

$(uj,uj) = S(0) = a(Q)- 1 dim R (H) 

for all w G O. 

On spheres, Condition (M) holds 10 for the spaces 

ft(0 (S"- 1 ) [respectively P^S™" 1 ), 
^(g"- 1 )] with S = [respectively S = C {1 \ S = R®]; the notation is that of 

Example 2.4. 

2.12. Proposition. Let (O, a), (F^^kyo, and (& k ^)k>o be as in Proposition 2.10. Choo- 
se I > and assume that Condition (M) holds for J 7 ^ and a function S®. 

If X is a tight geometrical 2l-design on O, the set of distances 

Dx = {c G ]0, oof | c = d(x, y) for some x, y G X, x ^ y} 

is contained in the set of zeros of the function . 

Proof. Notation being as in the proof of Proposition 2.10, we have 

V(x,x') = $W(x,x') = S^(d(x,x')) 

for all x,x' G X. Since ty(x,x') = when x ^ x' (see Example 2.3), the proposition 
follows. □ 

The last proposition of this section shows how the setting of reproducing kernels makes 
it straightforward to generalize a characterization which is well-known for spherical designs 
(Theorem 5.5 in [DeGoS-77], or Theorems 3.2 and 4.3 in [GoeSe-79]). 

10 Or rather the analogous condition for scalar products (lj | lj') rather than distances d(uj,uj'). 
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2.13. Proposition. Let (fi, a) be a finite measure space, let H be a finite dimensional 
Hilbert space of functions on O which is a subspace of L 2 (fi, cr), and let $ denote the 
corresponding reproducing kernel. Let X be a non-empty finite subset of Vt and let W be a 
weight function on X . 

(i) Wehavej: xyeX W(x)W(y)$(x,y) > 0. 

(ii) Equality holds in (i) if and only if Ylxex W(x)ip(x) — for all <p ETC. 

Proof. The inequality holds in (i) since $ is a kernel of positive type. 

Assume now equality. For any orthonormal basis (ej)j e j of H, Proposition 2.2.iv implies 
that 

W(x)W(y)*(x,y) = W{x)W{y)Y j eA^Av) 

x,yEX x,yEX jEJ 

= ^(^^(x) ej (x)) 2 = 0. 

jeJ xex 

Hence Ylxex W(x)ej(x) = for all j G J, and therefore Y^xex W(x)ip(x) = for all 
<pi "hi. 

The converse implication in (ii) is straightforward. □ 

2.14. Particular case. Consider two integers n>2 and k > 1, a non-empty finite subset 
X ofS n ~ 1 , a weight W : X — ► IR+ such that ^2 xeX W{x) = 1, and the polynomials 

RS^ defined in Example 2.4- Then 

Y W(x)W(y)Q ( - J \(x | y)) > and ^ W(x)W(y)R {j) ((x \ y)) > 1 

x,yEX x,yEX 

for all j > 0. Moreover, the following properties are equivalent: 

(i) (A, W) is a cubature formula of strength k on § n_1 ; 

(ii) equality holds in (*) for all j G {1, . . . , k}; 
(Hi) equality holds in (**) for j = k. 



3. The case of spheres 
Antipodal cubature formulas and spherical designs 

In this section, we assume that O = § n_1 is the unit sphere in R n , n > 2, that a is 
the rotation-invariant probability measure on § n_1 , and that the notation is as in Exam- 
ples 1.11 and 2.4. 

For tight (2/)-spherical designs, the set Dx of Proposition 2.12 is precisely known, and 
\Dx\ = I- Instead of distances and Dx, it is more convenient to use scalar products and 
the set 

A x = {c G [—1, 1[ | c = (x | y) for some x, y G A, x ^ y\ . 
The following proposition is Theorem 5.1 of [GoeSe-79]. 
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3.1. Proposition. If X is a tight spherical (21) -design on §> n 1 , then Ax coincides with 
the set of roots of the polynomial B> 1 ' defined in 2-4- 

Proof. Let denote the set of roots of R( l \ which is of order /. Since Ax C by 
the proof of Proposition 2.12, it is enough to show that a = \A X \ is not less than /. More 
generally, let us show that a > I for any spherical (2/)-design. 
Define for each x G X a polynomial function 7 X of degree a by 

lx ( u ) = Yl ^ 1 ^ - C for all uj G R n . 

ceA x 

Then r ) x (x') = 5 XjX ' for all x,x' G X, so that the family (•y x ) xeX is linearly independent. 
Thus 

diniR^CS"- 1 )) > \X\. 



bmce \X\ > dim M (J r W(§ n " 1 )) by Proposition 1.7, we have a > I. 

Ai) 

'R 



In particular, Ax = Z]9 for a tight spherical (2Z)-design. □ 



Propositions 1.7, 2.10, 2.12, 3.1, and Theorem 2.8 apply to cubature formulas and 
spherical designs of even strength. We expose now the analogous facts for odd strength. 

3.2. Proposition. (Compare with 1.7.) Let n > 2, I > 0, and P^S 71 " 1 ) c L 2 (S n ~ 1 ,a) 
be as in Example 1.11. Let moreover X = Y U (-Y) be a non-empty antipodal finite subset 
of the sphere and W : X — > be a symmetric weight. 

Then (X, W) is a cubature formula of strength 21 + 1 if and only if the restriction 
mapping 

p : V {1) — ► £ 2 (Y,2W) 
is an isometry. In particular, if (X, W) is a cubature formula of strength 21 + 1, then 

(«) \X\ >2dim K (P«(§- 1 )) = ^l^y 

Proof. Observe that ip(-uj) = (-l) j ip(u) for all j > 0, ip G P^S 71 " 1 ), and u G § n_1 , so 
that 

w(x)< fl (x)< f2 (x) = 2j2w(y)My)My) 

xex y eY 

for v?i,V2 eP'^S"" 1 ). 

The condition for p to be an isometry is 

2 V] W(y)(p 1 (y)(p 2 (y) = / Vi(w)V2(w)ti(T(w) 

^ /en— 1 

for all (pi, (p 2 G P^(§ n_1 ). As in the proof of Proposition 1.7, this can be written 



2 J2 W(yMy) = [ <p(u)da(u;) 
yev ^s- 1 
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or indeed 

(*) J2 W(xMx) = [ V (u)da(u) 

for all ip G V ( - 2l \S n ~ 1 ). Now (*) holds for all ip G p( 2 ' +1 ) (S™" 1 ), since both terms are zero 
in this case by symmetry reasons. Hence p is an isometry if and only if (*) holds for all 

v e jr(2I+l)(gn-l). n 



3.3. Remarks and definition, (i) Let (X, W) be an antipodal cubature formula on S n , 
and write X = Y U (-Y). Observe that, since a function on X can be canonically written 
as the sum of an even function and an odd function, we have a decomposition of l 2 (X 1 W) 
as the orthogonal direct sum of two copies of £ 2 (Y, 2W). It follows from Propositions 1.7 
and 3.3 that the restriction mapping 

jF«(§n-i) = p(O(s»-i) p(J-i)( S »-i) — ► £ 2 (X,W) = £ 2 (Y,2W)®£ 2 (Y,2W) 

is an isometry if and only if the restriction mapping pW(§ n_1 ) — >■ £ 2 (Y, 2W) is an iso- 
metry. 

(ii) For a spherical (21 + l)-design X in §> n_1 (not necessarily antipodal), it is known 
that \X\ > 2dim M (pW(§ n - 1 )); moreover, in case equality holds, then X is antipodal. See 
Theorem 5.12 in [DeGoS-77]; the proof uses the "linear programming method", which is 
more powerful than the "Fisher type method" used in our proof of Proposition 3.2. 

(iii) A cubature formula (X, W) of strength 21 + 1 on § n_1 is tight (see the definition 
in 1.11) if equality holds in Equation (ft) of Proposition 3.2. 

The following is Theorem 5.11 of [GoeSe-79]. 

3.4. Proposition. (Compare with 2.10, 2.12, and 3.1.) Let (X, W) be a cubature formula 
of strength 21 + 1 on §> n_1 which is tight. 

Then W is uniform, namely X is a tight spherical (21 + 1) -design. Moreover the set 

Bx = {c G ] — 1, 1[ | c = {x | y) for some x,y G X, x ^ ±y} 
coincides with the set of roots of the polynomial C"' defined in 2-4- 

Proof. Let zjp denote the set of roots of C^ l \ which is of order /. Since B x C by the 
argument used in the proof of Proposition 2.12, it is enough to show that b = \Bx\ is not 
less than /. More generally, let us show that b > I for any antipodal (21 + l)-design. 
Define for each x G X a polynomial function ^ x of degree b by 

7*M = JJ I ^ ~ c f or a ll 

ceB x ° 

As — Bx = Bx, we have 

n„ cR ((x \ uj) + c) , n cR ((x \ —u>) — c) , 
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for all uj E W 1 . Thus the restriction of % to S 71 " 1 is in ^(S™" 1 ). Decompose X as 
Y U (— Y) as in Item 1.12; since %(y') = 5 yjy > for all y, y' £ Y, we have 

dimM^^- 1 )) > \Y\ > dim K (P«(§ n - 1 )) 

and therefore b > I. 

In particular Bx = zj-) for a tight antipodal spherical (21 + l)-design. □ 

3.5. Remark. Given n and /, let (X, W) be a cubature formula of strength 21 + 1 on § n_1 
such that \X\ < \X'\ for any cubature formula of strength 21 + 1 on § n_1 . The weight W 
need not be uniform; see Example 3.7 for such a cubature formula of strength 7 on § 2 . 

3.6. Proposition. (Compare with 2.8.) For each I > 0, there exists a finite antipodal 
subset X = YU (-Y) o/S"" 1 such that 

dim^P^S"- 1 )) < \Y\ < dim K (P^(§ n - 1 )) -1 

and a symmetric weight W : X — > such that (X, W) is an antipodal cubature formula 
of strength 2/ + 1 on S n_1 . 

Proof: see Propositions 2.6 and 2.7. □ 

Similar estimates are known for other spaces, in particular for hypecubes in M n ; see 
page 366 of [DavRa-84]. 

3.7. Example: cubature formula of strength 5, 7, and 9 on § 2 . A number of explicit 
cubature formulas can be collected from the literature, either directly or indirectly. Many 
of them are reviewed in [Dicks-19] (see in particular pages 717-724) and [Rezni-92]. Let 
us first write down an identity of Lucas (1877): 

8(«i + u\ + u\) + ( M i + e 2«2 + e 3 u 3 ) 4 = I2(ul + uj + u 2 3 ) 2 . 

e 2 ,e 3 €{l,-l} 

Let Y C S 2 be the set of size 7 containing the three vectors ei,e 2 ,e 3 of the canonical 
orthonormal basis of M 3 and the four vectors 3 _1 / 2 (ei + €262 + £363), for €2, e 3 = 1, —1. 
Define a weight W : Y — > R* + by Wfa) = ^ and W(?,- 1 / 2 (e 1 +e 2 e 2 +e 3 e 3 )) = ^. Then 
Lucas' identity can be rewritten as 

V W(y)(y I w) 4 = -(u I u) 2 for all ueE 3 

yEY 

so that (y, W) is a cubature formula of size 7 for p( 4 )(§ 2 ) by the argument of Proposi- 
tion 1.13. If X = Y U (-Y) and if W is extended by symmetry, (X, \ W) is a cubature 
formula of size 14 and strength 5. 

Size \Y\ = 7 [respectively \X\ = 14] is not optimal for cubature formulas for V^(E> 2 ) 
[respectively JF( 5 )(§ 2 )]. Indeed, as already reported in Items 1.17 and 1.18, the 12 vertices 
of a regular icosahedron provide a tight spherical 5-design, and such a design is unique up 
to isometry [DeGoS-77, page 375] . The 6 vertices of a corresponding Y provide a design for 
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the space p( 4 )(§ 2 ). (Note that this is not a spherical 4-design, because ppj J2 y eY viv) * s not 
always equal to Jg 2 </?(w)<icr(c<;) for functions <p in the second summand of the decomposition 

JF(4) (§2) =p (4) (S2)0p(3)( S 2).) 

From [Rezni-92] and [HarSl-96] (see also Item 1.18), we collect the following facts. 

o There exists a pair (Y,W) with |Y| = 11 which is a cubature formula for p( 6 )(S> 2 ), 
and 11 is the minimal possible size (pages 133-135 of [Rezni-92]). Hence there exists a 
cubature formula (X, \W) for ^(S 2 ) of size \X\ = 22. On the other hand, a spherical 
design X C S 2 of size \X\ < 22 is of strength at most 5, and a spherical 7-design in S 2 
has size at least 24 (compare with the lower bound of Proposition 3.3 for 7-designs in § 2 , 
which is 20). 

o There exists a pair (Y, W) with \Y\ = 16 which is a cubature formula for V^ 8 \S 2 ), 
and 16 is the minimal possible size (pages 111 and 136 of [Rezni-92], referring to Finden, 
Sobolev, and McLaren). Hence there exists a cubature formula (X, \ W) for JF( 9 )(§ 2 ) of 
size \X\ = 32. On the other hand, a spherical design X C § 2 of size \X\ < 32 is of strength 
at most 7, and a spherical 9-design in S 2 has size at least 48 (compare with the lower 
bound of Proposition 3.2 for 9-designs in S 2 , which is 30). 

3.8. Example: cubature formulas of strengths 5, 7, 9 and 11 on S 3 . The following 
identities go back respectively to Liouville (1859), Kempner (1912), Hurwitz (1908), and 
J. Schur (1909): 

4 .4,2 

^(2w,) 4 + Yl {ui + t 2 u 2 +e 3 u 3 + e 4 u 4 f = 24f^w 2 J, 

i=l ei€{l,-l} M=l ' 

4 , 4 ,3 

J2(2u t ) 6 + 8j2(ui + euj) 6 + (u 1 + e 2 u 2 +e 3 u 3 + e 4 u 4 ) 6 = 120f^\ 2 J, 

i=l / ^€{1,-1} M=l ' 

4 

6^(2w,) 8 + 60^(w i + ew J ) 8 + ^2(2 Ul + ejUj + e k u k ) 8 

i=l I II 

+ 6 Y (ui+e 2 u 2 + e 3 u 3 +e 4 u 4 ) 8 = bOAofj^uf] , 

e<e{i,-i} m=i ' 

4 

9^(2^) 10 + 180 ^(w, + euj) 10 + J2(2 Ul + e jUj + e k u k y 



\z>Ui T tj Uj T tfcUfc llC 

i=l 7 77 



+ 9 ( Ul + e 2 w 2 + e 3 u 3 + e 4 u 4 ) 10 = 22680 (XX) ' 

e<e{i,-i} M=l ' 



Summations ^ 7 contain 12 terms ui + euj, with 1 < % < j < 4 and e = ±1. Summa- 
tions Yin contain 48 terms of the form 2ui + €jiij + e k u k , with k pairwise distinct in 
{1,2,3,4}, j < k, and €j,e k = ±1. These four identities provide cubature formulas on S 3 
of sizes and strengths 

24 and 5 (Liouville), 

48 and 7 (Kempner), 
144 and 9 (Hurwitz), 
144 and 11 (Schur). 
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Other examples appear in [Dicks-19] (pages 717-724). 

Observe that Liouville's identity provides a design for p( 4 )(§ 3 ) of size 12 and a spherical 
5-design of size 24. The latter is (up to homothety) a root system of type _D 4 ; in other 
coordinates, it can also be written as 



(L) 



(ui + euj) 4 = 6( 



U : 



l<i<j<4,e€{-l,l} 



i=l 



(apparently written down first by Lucas in 1876). This is not a design for ^^(S 3 ), 
otherwise there would exist a cubature formula of strength 7 on § 3 of size 24, but such 
cubature formulas cannot have size less than 40 by Proposition 3.2. 

Kempner's identity provides a spherical 7-design in S 3 of size 48. There is in [HarSl-94] 
a cubature formula of strength 7 in S 3 which is of size 46, conjecturally the optimal size. 

Schur's identity is not optimal for cubatures formulas of strength 11 on S 3 . Indeed, 
the 120 vertices of a regular polytope of Schlafli symbol 11 {3, 3, 5} provide an antipodal 
spherical 11-design (compare with the lower bound 112 of Proposition 3.3 for 11-designs 
in§ 3 ). 

Any orbit in § 3 for the natural action of the Coxeter group H 4 is a 11-design, and 
there exists a particular orbit of size 1440 which is indeed a 19-design. This construction 
is apparently due to Salihov (1975). See the end of Section 5 in [GoSe-81a], page 112 of 
[Rezni-92], as well as [HarPa-04]. 

3.9. A digression on Waring problem. Liouville used his identity (see 3.8) and La- 
grange's theorem on the representation of integers as sums of four squares to show the 
following claim: 

any positive integers is a sum of at most 53 biquadrates (= fourth powers). 
Here is a proof of the claim, using Lucas' form (L) of Liouville's identity. Since any 
positive number is (by Lagrange's theorem) of the form 6(Af + A| + A| + Af ) + r with 
A 7 !, . . . , A 4 G N and r G {0, 1, 2, 3, 4, 5}, it is enough to check that any number of the form 
6A 2 is a sum of 12 biquadrates. Using Lagrange's theorem again, A can be written as 
a sum n\ + n 2 . + n\ + n\ of four squares. If n G Z 4 denotes the vector of coordinates 
ni, n,2, n,3, 714, we have A = (n \ n) and 

6A 2 = 6(n I n) 2 = ^ (n^ + enj) 4 = sum of 12 fourth powers 

l<i<j<4,e€{-l,l} 

by (L), and this ends the proof. 

Today, we know how to make this bound sharp, reducing it from 53 to 19 [BDD-86a/b], 
or even to 16 if a, finite number of exceptions is allowed [Daven-39]. More precisely, there 
are exactly 96 numbers which are not sums of 16 biquadrates, and the maximum of them 
is 13 792 ; their list is shown in [DeHeL-00] . 

3.10. Lattice cubature formulas. Here is a variation on the construction of lattice 
designs described in Example 1.16. 



11 See [Coxet-73], page 153. 
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Consider a lattice AcK" which is even and unimodular. For an integer / > and a 
harmonic homogeneous polynomial function P G H ( ^ 2l ' > (W n ), the theta series is defined 
by 

e A , P (z) = ^p(a)^< a i a > = p(o)+ J2 p ( x )<i+ J2 p (x)q 2 + --- 

XeA x£A 2 x£A 4 

where q = e 2l7TZ , z G C with Im(z) > 0, and A 2 j = {x G A | (x \ x) = 2j}. It is a modular 
form for PSX(2, Z) of weight ^ + 21, and a parabolic modular form if moreover / > 0. 
In particular, we have 6a, p = whenever / > and ^ + 21 < 10 or ^ + 21 = 14 (see 
Corollary 3.3 in [Ebeli-94]). 

Suppose now that n = 8, so that A is a root lattice of type Eg. The previous conside- 
rations show that A 2 j is a spherical 7-design and a design for T^ 10 - 1 ^ 8 ) for any j > 1. 
Moreover, for P G 7i( 8 )(IR 8 ), the parabolic modular form Oa,p is necessarily a constant 
multiple of 

oo 

A 24 (z) = q 2 f[ (1 - Q 2m ) 24 = Q 2 - 24q 4 + 252q 6 - U72q 8 + ■■■ (q = e t7TZ ). 

m=l 

In particular, 24J2 xeA2 p ( x ) + E x eA 4 P <» = °- Observe that P(x) = 2 4 P(x/ v / 2) for 
x G A 2 and = 2 8 P(x/2) for a; G A 4 . It follows that 

x6A 2 \V^/ xeA4 

for all P G Ti^^M 8 ) and 1 < A; < 11. In other words, the first two shells of A provide a 
cubature formula in S 7 of strength 11, with underlying set -^A 2 U t?A 4 , of size 

|A 2 | + |A 4 | = 240 + 2160 = 2400, 

and with exactly two different weights, namely 1/1680 on -^A 2 and 1/2520 on |A 4 . This 
size compares favourably with the bounds of Proposition 3.6, which read here 

2dim M P (5) (R 8 ) = 1584 < N < 2dim M P (10) (M 8 ) - 2 = 38894. 

A similar lower bound shows that -^|A 2 U ^A 4 cannot enter a cubature formula of strength 
13 on § 7 . 

The previous construction indicates clearly a general procedure. 

3.11. On Gaussian designs. Let / be an integer, / > 0, and let a G R n . On the one 

hand, it is classical that 

Lj° 1 x)2 ' drT{x) = n(n + 2 ( ) 2 '>"a-2) ( " 1 °>' 

(notation: (2/— 1)!! = Y\ l j = i(2l~2j+1)). Indeed, the left hand side is clearly a homogeneous 
polynomial of degree 21 in a which is invariant by the orthogonal group 0(n), and therefore 
a constant multiple of (a \ a) 1 . If we apply both sides the Laplacian with respect to a, 
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we find a recurrence relation for the constants, and this provides the result. On the other 
hand, we have 

1 f ,.,„v2« -INI 2 ,™ _ (21-1)", , w 



(a I xY l e-" x " dx = v ' ' (a I a) 1 
(calculus, using integration by parts in case n = 1). It follows that 

(*) / ^(x)dcr(x) = — — — % — — — ^ [ <f(x)e~^ 2 dx 

V ' y s «-i re(ra + 2) • • • (ra + 2Z - 2) W 2 7 

for all </? G p( 2 ')(M n ), / > 0. [See the footnote in Proposition 1.13.] Observe that, for 
tp G "P( 2Z+1 )(IR n ), both integrals in (*) vanish for symmetry reasons. 

Let now X be a non-empty finite subset of of size N. If X is a spherical t-design 
with t G {1,2,3}, the set 

is a Gaussian t-design. Indeed, for (p G V^CMJ 1 ), we have 



— ip(x) = (fix) = — ip(x)da(x) = — 777 / ¥>(x)t 



|2 



(we leave to the reader the verifications with <p G "P (1) (R n ) and <^ G "P (3) (R n )). 
If X is now a spherical 5-design, set 



Pi = sJ\{n-V2n) \{p\ + pl) = \ 

, so that " . 

n 1/4 4n n(n + 2) 

p 2 = y 2 (n + ^ 2 (pl + p2) = — i — • 

Then 

Xgg = Pi^ U P2-^ (disjoint union) 
is a Gaussian 5-design of size 2iV. Indeed, for (p G p( 4 )(R n ), we have 

(where the last equality is a consequence of the values chosen for p\ and p?). Similarly 
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for all ip G V^'(M. n ). We leave to the reader the verifications with (p a homogeneous 
polynomial of odd degree. 

We end the present chapter with a proposition and a construction suggesting one more 
connection with another subject. 

For an integer N > 1 and a real number p > 1, we denote by £ P (N) the classical Banach 
space of dimension N, with underlying space M, N and norm ||x|| , where 



( x j)i<j<N — \l 



\ x j 



l<j<N 



A linear mapping from a real vector space to £ P (N) is said to be degenerate if its image 
is in a hyperplane of equation Xj = for some j G {1, . . . , N}. The next proposition is 
from [LyuVa-93], and en is the constant of our Proposition 1.13. 

3.12. Proposition. Consider integers n > 2, I > 2, and N > 1. 

(i) To any cubature formula (X, W) of size N for V^ 2l) (S 71 - 1 ) corresponds a non-dege- 
nerate isometric embedding Jx,w '■ £ 2 (n) — > £ 2l (N). 

(ii) Conversely, to any non- degenerate isometric embedding J : £ 2 {n) — > £ 2l (N) corres- 
ponds a cubature formula of size N in V^i^- 1 ). 

Proof. Let (X, W) be as in (i) and let (x^\ . . . , x^) be an enumeration of the points in 
X. The linear mapping Jx,w '■ ^ n — ► ^ N defined by its coordinates 

^ M (x {k) \u), l<h<N, 

is an isometry £ 2 {n) — ► £ 2l (N) by the argument of Proposition 1.13. It is non-degenerate 
since X 1 - = {0}. 

Let J : £ 2 (n) — ► £ 2l (N) be a non-degenerate isometry as in (ii). For each k G 
{1, . . . , N}, let y^ be the unique vector in £ 2 {n) such that J{u)k = (y^ \ u) for all 
u G £ 2 (n); observe that y( k > ^ 0, by the non- degeneracy condition on N. Set x^ k ' = 
y {k) /\\y (k) \\ 2 and = C2l \\y^\\l l . Then X = {x^\ . . . , i^} is a A-subset of S™" 1 , 

the define a weight X — > R+, and (A, W) is a cubature formula for £>( 2 (§«•-!) by 
Proposition 1.13. □ 

4. Markov operators 

Given a complex Hilbert space H, we denote by C(H) the C*-algebra of all bounded 
linear operators on 7i and by U(7i) its unitary group, consisting of all unitary operators 
on 7i. Given a group T, we denote by C[r] the group algebra of all functions V — > C 
of finite support, with multiplication the convolution; we consider V as a subset of C[r], 
by identifying 7 G V with the function of value 1 on 7 and elsewehere. A unitary 
representation of a group V in a Hilbert space Ti is a group homomorphism n : V — > U (H) ; 
it extends to a morphism of algebra C[r] — > >C(7Y), denoted by n again, and defined 
by tt(/) = S 7 er fil^il)- Recall that the norm of a unitary operator is 1, so that 

IK(/)II < E 76 r 1/(7)1- Reca11 also that ' if /(7" 1 ) = /(t) for a11 7 ^ r, the operator tt(/) 
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is self-adjoint; in particular, the spectrum of 7r(/) is then a closed subset of the real interval 

l-\W)h\W)\\}- 

Up to minor terminological changes, this holds for a real Hilbert space 7Y, its ortho- 
gonal group 0(TC), orthogonal representations V — > 0(H), and corresponding morphisms 



4.1. Definition. Let V be a group, S a finite subset of V, and W : S — > a weight 
function; set M$,w = X] S 6S^ /r ( s )' s e K[r]. The Markov operator associated to the 
pair (£, W) and an orthogonal (or unitary) representation it : V — > is the operator 

tv(M SjW ) = J2w(s)tt(s) 

ses 

on 7i. 

Markov operators (often with constant weights) play an important role in the study of 
random walks on Cayley graphs of groups [Keste-59], and more generally in connection 
with unitary representations of groups. See for example [BarGr-00], [BeVaZ-97], [GaJaS- 
99], [GriZu-01], [GriZu-02], [GriZu], [HaRVl-93], and [HaRV2-93]. 

Consider a sigma-finite measure space (O, a) and a group G acting on O by measurable 
transformations preserving the measure class of a. The corresponding unitary representa- 
tion p of G on L 2 (0, a) is defined by 



(p(g)ip)(uj) = yj^- <p(g 1 u) 



for g G G, ip G L 2 (0, a), and u G O, where denotes the appropriate Radon-Nikodym 
derivative. Given a mapping s : O — > G, u i — > s w , a finite subset X of O, and a weight 
function W : X — ► M*^_, we have an element Mx,w = Ylxex W{x)s x in the real group 
algebra of the subgroup of G generated by s(X). Given moreover a representation tt of G 
in a Hilbert space 7i n , we have a Markov operator 

xex 

The particular case studied in [Pache-04] is that of the orthogonal group G = 0(n) act- 
ing on the unit sphere O = S n_1 , with one of the finite-dimensional spaces introduced 
in Example 1.11 above. In this case, the image of a point x G § n_1 by the mapping s is 
the orthogonal reflection of M n which fixes the hyperplane orthogonal to x. 

4.2. Proposition. Let n > 2, I > be integers and let denote the natural represen- 
tation of the group O(n) in the space Ti,^(E> n ~ 1 ) of harmonic homogeneous polynomials 
of degree I. For each x G § n_1 , let s x G 0(n) denote the reflection of M 71 that fixes the 
hyperplane x of~R n . 

If (X, W) is a cubature formula of strength 21 on then 

n^(M xw ) = idffl 
v *,wj 2l + n-2 
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where i 

d (0 

denotes the identity operator on H^(S n 1 ). 

Proof, repeated from [Pache-04] (see also [HarVe-01]). For x G M n , we define the operator 
s x on M n by s x (u) = (x \ x)u — 2(x \ u)x. Note that s x is selfadjoint, that §\ x = X 2 s x for 
AgI, and that s x = s x if x G S n_1 . 

For k > 0, define the selfadjoint operator 



M {k) = [ ^ k \ SuJ )da{uj) 



on 7i^ fc )(M n ). For g G O(n), we have s g ^ = gs w g 1 . Since the measure a is 0(n)- 
invariant and the representation is irreducible, it follows from Schur's lemma that 
is a homothety. By a simple computation, we obtain the trace of ^^(s^), which is 
independent of u;, and therefore the trace of M .It follows that 

(1) M (fc) = n ~ 2 id<*> 

for all k > 0. 

For v?i,v?2 G ft( fc )(M n ) and u G R n , set 

^ (fc) (</?!, 992) («) = (<£i I ¥2 ° s«) = / ip 1 {u)ip 2 {s u {u))do-{u). 

It is easy to check that the function ty( k \ipi, 922) is polynomial, homogeneous of degree 2k, 
and depends symmetrically on <pi,<f2- Thus, we have a linear mapping 

: Sym 2 (ft (fc) (R n )) — > P (2fc) (R n ). 

Moreover ^.^Wf^.^jH^H = ^ | M ( % 2 )). 

Consider now a non-empty finite subset X of § n_1 , a weight : X — > , and an 
integer / > 0. Then 

(2) > W(x)ip(x) = / (p(u)da(u) for all y> in the image of ' 



^ W{x){ Vl I v?2 o s K ) = U x M ( % 2 )) for all (fx, ^ 2 G 7i (z) (R n ), 



if and only if 

— (/) 

xex 

namely, by the equality tt"> (Mx,w) ( P2 = J2 xe x W(x)(p2 s x defining tt">(Mx,w) and 
by (1), if and only if 

(3) *®{Mx, W ) = M W = ^^id«. 

If (X, W) is a cubature formula of strength 21, then (2) holds for all ip G V^ 2l \M. n ), and 
a fortiori for all ip in the image of ^ l >, so that the proposition is proved. □ 



It is remarkable that a converse holds for n > 2. 
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4.3. Proposition. Let n > 3, / > be integers and let tt^ be as in the previous propo- 
sition. Let X be a non-empty antipodal finite subset o/S n_1 and let W : X — > be a 
symmetric weight. 

If 7fO (Mx,w) is a homothety, then (X, W) is a cubature formula of strength 21 + 1. 
Proof. See [Pache-04]. □ 

More generally, consider a Riemannian symmetric pair (G, H) where G is a compact 
Lie group acting on O = G/H, and the G-invariant probability measure a on O. Let 
s = s symrn : O — > G be the mapping which associates to a point x E G/H the symmetry 
of O fixing x. (For spheres, observe that s symm (x) is minus the reflection s x fixing x^ .) 
There is an orthogonal decomposition 

L 2 (0,a) = 07 A 
AeA 

in irreducible G-spaces. For a finite subset X of Vt and a weight function W on X, there is 
an analogue of Proposition 4.2 concerning spaces V x for which [X, W) is a cubature formula 
and for which the Markov operator -k^ x \Mx,w) £ End(V^ A )) is a constant multiple of the 
identity [Pache, in preparation]. For designs in G/H a Grassmannian, see [BaCoN-02] and 
[BaBaC-04]. 

Consider a non-empty finite subset S of 0(n) which is symmetric (S -1 = S) and denote 
by Ts the subgroup of 0(n) it generates. 

From now on, we assume that W(s) = | *S* | 1 for all s £ S 
and we set M s = {Sf 1 J^seS s G R [ T s}- 
Propositions 4.2 and 4.3 relate the spectra of Markov operators Tx^ k \Ms) for some small 
values of k to the cubature properties of (X, W). It happens that the spectra for k — > oo 
of the operators n^(Ms) are also important, as shown by Lubotzky, Phillips, and Sarnak 
[LuPhS-86, 87] (see also [ColdV-88]). We formulate the weak Observation 4.4 before the 
stronger Proposition 4.5. 

Let 7To denote the natural representation of 0(n) in the space L^S 71-1 , a) of L 2 -functions 
of zero average. The orthogonal sum Lo(S> n_1 , a) = ®^ =l 1iS k ^ (§ n_1 ) provides an or- 
thogonal decomposition n = ®^ =1 k^ into irreducible subrepresentations (see, e.g., 
[SteWe-71]), so that 

4.4. Observation. For any finite symmetric subset S ofO(n), we have 



1 > ho(Ms)\\ > 




Proof. (Compare with Proposition 2.3.2 in [Sarna-90].) For all s G S 1 , we have ||7To(s)|| = 1, 
since the representation n is orthogonal. The upper bound follows. 

Choose io G S n_1 such that s(u) ^ t(u) for all s,t e S U {id}, s ^ t. There exists a 
neighbourhood U of u in § n_1 such that s(U) n t(U) = for all s,t e S U {id}, s ^ t, 
and a function cp e Lo(S n_1 ,(7) of norm 1 supported in U. Since ||7To(Ms)</?|| 2 = \S\ , we 
have |MM S )|| > |S|~ 1/2 . □ 



7Tq(Ms)\\ = SUp 
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4.5. Proposition. Let S be a symmetric finite subset of 0(3). Then 



|7To(Ms)|| > limsup 

k^oo 



> 



The proofs of [LuPhS-86] and [ColdV-88] are written up for the case of a subset S of 
SO(n). The present generalization to 0(n) is rather straightforward. We isolate part of 
the proof in the following lemma 

4.6. Lemma, (i) Let g G SO(3) be a rotation of angle 9 g G [0, n]. Then 

f sin[(2fc + l)fl g /2] 
trace = sin^/2] * J 

[ 2k + 1 if 6 g = 0. 

(ii) Let g G 0(3) have eigenvalues exp(±i9 g ) and —1. Then 

( cos[(2fc + ljgg/gj 
trace (*<*>(</)) = i cos[^/2] 7 s ^ 

( (-l) fc i/ ^ = 7T. 

Proof, (i) The eigenvalues of g G 5'0(3) are exp(±i6> 9 ) and 1. The space p( fc )(S 2 ) has a 
linear basis of eigenvectors of the transformation induced by g of the form 

{x a y b z k - a - b | a, 6 > and a + b < k) . 

For j G {0, . . . , k}, the trace of the linear endomorphism defined by g on the linear span 
of {x a y b z k - a ~ b | a, b > and a + b = j} is 



sin[6> 3 ] 

so that the trace of the linear endomorphism defined by g on "P( fc )(§ 2 ) is Yl > j=o (T j(^g)- 
Since P( fc )(§ 2 ) = e^o^-^o^" 2 ^ 2 ), we have trace (*-(%)) = a fc (^) + a k _ x (9 g ) 
and the formula of (i) follows. 

(ii) Similarly, for g G 0(3), (7 G" £0(3), the trace of the linear endomorphism defined 
by g on P«(§ 2 ) is £j =0 (-l) fc_i <7j ( 9 g) so that 

trace (tt^)) = a k (6 g ) - a k -i(O g ). 
(The only part of the proof which is not explicitely in [LuPhS-86] is (ii).) □ 

Proof of Proposition 4-5. Step one. Recall that T$ is the subgroup of 0(3) generated 
by S. For each integer N > 0, let Wn denote the number of words 12 in letters of 5, of 



12 Words which need not be reduced in any sense. 
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length N. To each word w G Wn corresponds naturally an orthogonal transformation 
g = g(w) eF s C 0(3). We have (M S ) N = \S\~ N E we w N dW and therefore 



\S\» trace UW((M a )«j) = V ^±#^ + V ^ + ^ 
11 V vv V ^ sin[0„/2l ^ cosf#„/2l 

by the previous lemma. 



sin[0„/2] " cos[#„/2] 

9W6SO(3) g (r»)^SO(3) 



Let mjv denote the quotient by \S\ N of the number of words w G W^r such that g(w) = 



1 G 0(3). We have 

for each N > 0. (This is Theorem 1.1 in [LuPhS-86]; observe however the change in 
notation: tjin there is \S\ N times what mjv is here and in [Keste-59].) 

Step two. Consider the Cayley graph 13 Cay s of Ts with respect to S, the left regular 
representation As of Ts, the corresponding Markov operator Xs(Ms), and its spectral 
measure /is- Kesten has shown that 

m N = J R t N d/j s (t) for all N > 0, 

||As(M s )|| =limsupjv =00 

\\Xs(Ms)\\ > 2y/\S\ — 1/ \ S\ with equality if and only if Cay s is a tree. 
See [Keste-59]. 

Step three. For each k > 0, let jis,k denote the quotient by dim^T^^S 2 ) = 2k + 1 of 
the spectral measure of yr( fc ) (M s ); this can be written 



^ 2k+l 



where 5(A) denotes a Dirac measure of support A and where (Afc )J ) 1<J<2A ,_ | _ 1 are the eigen- 
values of the endomorphism n^ k '(Ms) of ^^(S 2 ). Let N > 0; from the definition of ns,k, 
we have 

^J— trace ^ k \(M s ) N )) = jfVdps,* ; 
from the two previous steps, we have 

lim — L- trace (^ k \(M s ) N )) = m N = [ t N df, s ; 

it follows that the sequence of measures (/J,s,k) k>1 (all of mass at most 1 since ||7r( fc )(Ms)|| < 
1 and ||As(Ms)|| < 1) converges vaguely to us, and in particular that 



lim sup 

fc^oo 



7T (fc) (Ms) 



2y^^T 

\S\ 



□ 



13 Two elements 71, 72 G Tg viewed as vertices of Cay s are joined by an edge if 7 1 1 72 G S. We assume 
that 1 ^ 5, so that the Cayley graph is simple (= without loops), of degree 
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4.7. Theorem (Lubotzky-Phillips-Sarnak). For each prime p, there exists a subset 
So of SO (3) of p + 1 elements such that, if S = Sq U (Sq) -1 , 



n 



(k) 



(M s ) 



< 





s 


- 1 




s 





P+i 



for all k > 1. 



4.8. Open problem. It is a natural problem, for each / > 0, to look for a finite set 
S C 0(3) of reflections (depending on /) such that the following conditions are fulfilled: 



(*) 



TX {k \M S ) 



2k + 1 



id (fc) for k G {!,...,/} 



(see Propositions 4.2 and 4.3), 



< 



2^/\S\^l 



for k large enough 



(see Proposition 4.5 and Theorem 4.7), and \S\ as small as possible. 

Observe that (*) is a way to write that all eigenvalues of tc^ (Ms) are equal to l/(2/c+l), 
for k G {1, ... , /}, and that (**) is a bound on the eigenvalues of (Ms), for large values 
of k. 



5. Reflection groups 



Let V be a Euclidean space. The reflection s x G 0(V) associated to x G V, x ^ 0, is 
defined as above by 

s x (y) = V ~ 2j-\^-x for all y G V. 

Observe that s x > = s x if and only if Rx' = Rx. In this chapter, a reflection group is a 
subgroup of OiV) generated by a finite set of reflections; it need not act properly on V, 
contrarily to what is assumed most often in [Bourb-68]. Thus, any finite subset X of 
V \ {0} defines a reflection group Fx C O(V) generated by (s x ) x6X . 

In particular, let A be a lattice in V. We assume that A is integral, so that A is the 
disjoint union of the origin and of its non-empty shells A m , m > 1, as defined in Item 1.15. 
Let rA, m denote the reflection group generated by {s x | x G A m }. For example, if A is 
a root lattice (namely an even lattice A generated by A2), then Ta,2 is a finite group, 
and indeed a direct product of Coxeter groups of types A n (n > 1), D n (n > 4), and E n 
(n = 6, 7, 8); see for example [Ebeli-94]. In most cases however, TA )m is an infinite group. 

Let us specialize to m a power of 2. The scalar product on V defines a symmetric 
Z-bilinear form on A which extends to a symmetric Z [1/2] -bilinear form 

: (A (8) Z[l/2]) x (A <g> Z[l/2]) — ► Z[l/2]. 

We denote by 0(A, Z[l/2]) the orthogonal group of Observe that r A , m C 0(A, Z[l/2]). 

In case of the cubical lattice Z n C M n , we write 0(n, Z[l/2]) rather than 0(Z n , Z[l/2]). 
Several lattices A define the same group 0(n, Z[l/2]); a sufficient condition for the group 
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0(A,Z[l/2]) to be isomorhic to 0(n,Z[l/2]) is that 2 fc Z n C A C 2~ l Z n for some fc,Z > 0. 
In particular, if A C R 8 is a root lattice of type E 8 , then 0(A, Z[l/2]) = 0(8, Z[l/2]), since 
A and Z 8 have a common sublattice of index 2 (which is a root lattice of type Dg). Here 
are some properties of the group 0(n, Z[l/2]): 

(0) it is generated by reflections for any n; 

(1) it is infinite if and only if n > 5; 

(ii) it is naturally 14 a discrete cocompact subgroups of the 2-adic group 0(n, Q2); 

(iii) it is finitely presented; 

(iv) it is virtually torsion free, and more precisely Ker (0(n, Z[l/2]) — >■ 0(n, Z/3)) is 
torsion free; 

(v) it has finite-type homology — more precisely, for any Noetherian ring A and any 
i > 0, the A-module Hi(0(n 7 Z[l/2]), A) is finitely generated; 

(vi) it is of virtual cohomological dimension |(n — inf m6 ^ \n — 8m\), and the group 
Hi(0(n, Z[l/2]), Q) is zero in dimensions neither and nor maximal; 

(vii) the natural inclusion of 0(n, Z[l/2]) into the compact orthogonal group O(n) has 
dense image if and only if the group 0(n, Z[l/2]) is infinite. 

Several of these properties are straightforward consequences of (ii). Property (vii) follows 
from a strong approximation theorem due to Kneser [Knese-66] . For much more on these 
groups, see [Colli]. 

The next proposition shows a remarkable generating set of 0(8, Z[l/2]); it is an unpub- 
lished result of B. Venkov. 

5.1. Proposition (B. Venkov). Let L be a root lattice of type Eg. Then the group 
0(8, Z[l/2]) is generated by Tl,a and a reflection with respect to one root, and also by the 
finite Weyl group Vl,2 of type Eg and a reflection with respect to one element of L4. 

Proof. Let C be the set of even unimodular lattices M mV which are such that M(~)L is of 
index a power of 2 in both L and M (the same power since L and M are both unimodular). 
Then £ is a metric space for the distance 5 defined by 

5(M, N) = log 2 ([M : M n N]) = log 2 ([N : M f] N]). 

It is a theorem of Kneser [Knese-57, page 242] that £ is connected by steps of length 
one; more precisely, given M,N e C with d = 8(M,N), there exists a sequence M = 
M(°), MW, . . . , M^) = iV in C such that S(M^~ 1 \ M^) = 1 for all j G {1, . . . , d}. 

Set G = 0(8, Z[l/2]); this group operates naturally on C. Here is the synopsis of the 
proof: for g E G, set M = g(L) G C and d = 5(L, M); we show by induction on d that g 
is a product of appropriate reflections. More precisely, for = L, . . . , = M 

as above, we check that = s x (M( j ~^) for some x G M^~^ with (x | x) = 4. 

Step one: (T Ly2 , V L ^) = (s r , r L;4 ) = (T Lj2 , s x ) for any r G L 2 and x G L 4 . 

(For a subset S C 0(V), we denote by (S) the subgroup of 0(V) generated by S.) 



14 More generally, if P is a finite set of distinct primes and if iV denotes their product, 0(n, Z[l/N]) 
is naturally a discrete cocompact subgroup of the product of 0(n) with IlpeP ^( n ' Qp)> an d therefore 
also a discrete cocompact subgroup of Ilpep 0(n, Q p ). It follows that 0(n, Z[l/N]) is finitely presented 
[Serre-71]. 
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It is known that the group of all automorphisms of L coincides with the group ^l,2 
generated by the symmetries with respect to the roots. For g G Aut(L) and x G L, x 7^ 0, 
we have 

(*) s 3 (x) = gsxg' 1 - 

As Aut(L) acts transitively on L 4 , we have (T Ly2 , ^l,4) = (Tl,2, s x ) for any x G L 4 . (It is 
known that, more generally, Aut(L) acts transitively on each of L2, L 4 , Lq, the complement 
of 2L2 in L§, L10, and L12; see page 122 of [ConSl-99].) 

Relations (*) show that T L2 H ^l,4 is a normal subgroup of • Let ei, . . . , es be an 
orthonormal basis of V such that r = e\ + e-i and r' = e\ — e<i are in L2 (see page 268 of 
[Bourb-68]); then x = 2e\ and x' = 2e2 are in L 4 . A straightforward computation shows 
that idy 7^ s r s r > = s x s x > G r L)2 nr ii4 . As T L ^ is almost simple (by Exercise 2 of § VI.4 in 
[Bourb-68]), it follows that T^ 2 H 4 is a subgroup of index 15 at most 2 in so that 
Tl,2 is generated by r L)2 n T LA and s r for any r e L 2 . A fortiori (1^,2, r Li4 ) = (s r , T Lj4 ) 
for any r G £2- 

S^ep ^0; a reminder on neighbours. 

For this step, V can be a Euclidean space of any dimension n = (mod 8) and M any 
even unimodular lattice in V. We view M = M/2M as a vector space of dimension n 
over the prime field F2. There is a nondegenerate quadratic form q : M — > ¥2 defined by 
q(z) = \{z I z) (mod 2) for z G M representing ~z G M. For each z G M such that 2; ^ 2M 
and (z I z) = (mod 4), set 

M z = {m G M I (m | z) = (mod 2)} and M 2 = M z U (jL + M z ). 

Then M z is an integral unimodular lattice in V, and M n M z is of index 2 in both M and 
M z . Moreover: 

(i) M z is even if (z | z) = (mod 8) and odd if (z | z) = 4 (mod 8); 

(ii) for two ^, 2/ in M, not in 2M, with (z | z) = (mod 4) and (z' \ z') = (mod 4), 
M 2 = M 2 ' if and only if z' - z G 2M and (z | z) = (z' | z') (mod 8); 

(iii) any integral unimodular lattice M' in V such that M D M' is of index 2 in both M 
and M' appears as one of the lattices M z (z G M, 2; ^ 2M, (2; | z) = (mod 4)). 

(For the analogous facts concerning an odd lattice M, see [Venkov-79] ) . 

Step three: short representatives in M of non-zero isotropic classes in M. 

Let V be again of dimension 8, so that two even unimodular lattices in V are always 
isomorphic, and let M be such a lattice. The 2 8 elements of M splits as 

(a) the origin, 

(b) 120 elements represented by pairs ±r of roots, which are nonisotropic for q, 

(c) 135 nonzero isotropic elements. 

Let ip : M4 — > M be the restriction to M 4 of the canonical projection. We claim that 
each fiber of ip has at most 16 elements. Indeed, for x\, xi G M 4 such that x\ 7^ X2 and 
xi = x~2, there exists m G M such that X2 — x\ = 2m 7^ and, upon changing X2 to — X2 
if necessary, {x\ \ X2) > 0. Then 

(i) (x 2 - x\ I X2 — x\) = 4 + 4 - 2{xi I X2) < 8, 

(ii) (x2 — x\ I X2 — xi) = 4(m | m) > 8, 



15 We do not know if T^ ^ PI 4 is the whole of ^l,2 or if it is the subgroup of elements which are 
products of even numbers of reflections. 
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so that (x2 — %i | %i — X2) = 8 and x\ _L x<i- The claim follows. 

Since M 4 has exactly 2160 = 135 x 16 elements and since the image of ip is a priori 
inside the set of 135 nonzero isotropic elements of M, it follows that ip is onto this set. 

Step four: Let M, M' cV~K 8 be two even unimodular lattices such that M fl M' is of 
index 2 in both M and M' . Then there exists x G M4 such that M' = s x (M). 

Choose a nonzero isotropic element ~z G M represented by an element x G M4 (see Step 
three). Since M is unimodular, there exists a root r G M2 such that (r | x) = 1. The 
element z = x — 2r is also a representative of z, and 

(z\ z) = (x I x) - 4(r ] x) + 4(r | r) =4-4 + 8 = 8. 

Thus M z is the even neighbour M' of M (see Step two). 

To conclude the proof of Step four, we have to check that s x (M) = M z . For m G M, 
we have 

s x (m) = m — —(z + 2r | m)(z + 2r) 

= m — -(z I m)z — (z I m)r — (r \ m)z — 2(r | m)r 

and 

(z I s x (m)) = (z \ m) — (z \ m) (z \ r) (mod 2). 

If (z I m) is even, namely if m G M z , these formulas show that s x (m) G M z , and con- 
sequently that s x (m) G M z . If (z \ m) is odd, namely if m G M \ M z , they show that 
s x (m) — \z G M 2 , and we have again s x (m) G M 2 . Thus s x (A^) C M z . Since s x (M) and 
M z are both unimodular, this shows that s x (M) = M z . 
[Conversely, for any x G M 4 , we have 8{M, s x (M)) = 1.] 

End of proof of Proposition 5.1. 

Consider, as in the beginning of this proof, an element g G (9(8, Z[l/2]), the lattice 
M = g(L) G C at distance d from L, and a sequence <d such that = L, 

= M, and M^ -1 ) n M^') of index 2 in both M^" 1 ) and M {j \ The previous steps 
show that there exists for each j G an element G (M^ J_1 ^ 4 such that 

s Xj (M( J_1 )) = M"( J ). Set <7j = Since g~ 1 gd is in Aut(L), and in particular in 

(9(8, Z[l/2]), by Step one, it is enough to show that g d G 0(8,Z[l/2]). We claim that 
gj G 0(8, Z[l/2]) for all j G {1, . . . , <i}, and we prove the claim by induction on j. 

The claim is clear for j = 1. Assume it holds for some value oi j G {1, . . . , d — 1}, and 
let us check it for j + 1. Let y 7 - G L be the element such that gjiyj) = xj G Then 
(7 J+ i = s Xj gj = gjS Vj g~ 1 gj is indeed in G since G G and s^. G r^^. This shows that 
the claim holds for j + 1. □ 

5.2. Corollary. Let L be a root lattice of type E$. Then the subgroup SO(8, Z[l/2]) 0/ 
elements of determinant +1 in 0(8,Z[l/2]) is generated by the products s x s x i , where x,x' 
are in the shell L4 of L. 

5.3. A set of 8 generators for 0(8,Z[l/2]). Let (ei,...,e 8 ) denote the canonical 
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ortho normal basis in M 8 . As in [Bourb-68], set 

«i = -^{ei + eg) - -(e 2 H h e 7 ) 

"2 = ei + €2 

«i = Cj-i - Cj-2 (j = 3, ...,8) 

so that {«i, . . . , a 8 } is a basis of a root system of type E$, and therefore also a basis of a root 
lattice L C M 8 of type i? 8 - Let si, . . . , s 8 denote the reflections associated to cti, . . . , a 8 . 
Then Aut(L) has a Coxeter presentation with generators s±, . . . , s 8 and relations of the 
familiar form (s/sj )'"■■■• = 1. 

Consider the vector 2e2 G L4 and the corresponding reflection §2 : x 1 — > x — (x | 62)62- 
A simple computation shows that the conjugation by §2 exchanges S2 with S3 and leaves 
sj invariant for j = 1,4,5,6,7,8. It follows from Proposition 5.1 that 0(8, Z[l/2]) has a 
generating set obtained from that of Aut(L) by replacing S2 by §2- Moreover, the order 
m2,j of §2Sj in 0(8, Z[l/2]) is equal to 

00 for j ' = 1 
4 for j = 3, 4 
3 for j = 5, 6, 7, 8. 

Thus, for the generators of (9(8,Z[l/2]) described here, the orders rhij of the products of 
two generators coincide with the corresponding orders rriij for the Coxeter generators of 
Aut(L), this for all but three pairs of indices, namely for all but (2, 1), (2, 3), and (2, 4). 
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